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The Classical Theories of Radiation Reaction*) 


THomas ERBER 
Illinois Institute of Technology, Chicago 


This article presents a critical review of the classical theories of radiation reaction. The 
renormalized Dirac theory is discussed in particular detail. It is shown that all regularization 
prescriptions for this theory may be derived from a general dynamical correspondence prin- 
ciple. It is also shown that the physical content of the theory is severely limited by a general 
radiation condition. Various non-local and extended electron theories are classified with the 
help of the Herglotz-Wildermuth runaway‘ theorem. We discuss the existence of electro- 
dynamic collective modes; in particular the radiationless and self-excited states of extended 
charge structures. A possible connection with the existence of the u-meson is pointed out. 
In the concluding section various problems raised by the phenomenological radiation reaction 
theory of Ginzburg and Eidman are considered. 
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1. Introduction 


The fundamental premise of all theories of radiation reaction is that there are some 
dynamic consequences inherently associated with every elementary act of radi- 
ation. If this concept is incorporated into the structure of ordinary electrody- 
namics then the resulting theory can be conveniently summarized in terms of 
the following block diagram: 


Steering and Particle 


Radiation 


accelerating fields dynamics 


The links Q) and () here represent the ordinary Lorentz force and the Liénard- 
Wiechert potentials. The effects of radiation reaction may in principle then be 
taken into account through the insertion of the feedback“ loop indicated by @) . 
In the overwhelming majority of practical cases the coupling to the radiation 
field may be considered sufficiently weak so that any iterative corrections through 
this feedback loop can be safely disregarded. There are however certain excep- 
tional circumstances, for instance the self-excited states characteristic of electro- 
dynamic collective modes (see Section 6-G), in which the cycling through the 


closed loop QE is indeed the decisive element. 


Quite apart from any direct experimental motivation, the idea of radiation 
reaction also enters into conventional field theories through the logical necessity 
of balancing the energy-momentum conservation laws. The radiation reaction 
theories considered in Section 2 are in fact all constructed on this basis. The 
essential point here is that the energy and momentum losses computed directly 
from the Liénard-Wiechert potentials (i. e. @ above) may be exactly compen- 
sated through the addition of a suitable reactive force (Freact) to the ordinary 
Lorentz force law. Such a procedure has the attraction of being almost completely 
model-independent and also is easily adapted to various covariant generali- 
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zations. However there always remains an essential gap in this type of argument 
(Section 6-B) since the conservation laws of course do not uniquely specify any 
particular equation of motion. 

A completely different line of approach, which avoids all these ambiguities, is 
explored in the various non-local and extended electron theories studied in Sec- 
tion 6. The basic point of departure in these cases is the observation that Newton’s 
third law of action and reaction does not strictly apply to retarded action-at-a- 
distance theories. If for example the electron is visualized as being a sort of 
extended electromagnetic sponge, and if the validity of the Maxwell-Lorentz 
theory is extrapolated down to infinitesimal elements in its interior, then it is 
indeed possible to show that the net uncompensated intra-particle forces can give 
rise to ’’self-dragging“ or radiation reaction effects. 

Various modifications of these two major types of radiation reaction theories are 
discussed in several other sections of the paper. The principal conclusions of the 
present review are summarized in the paragraphs (i) — (iv) listed below. 


(i) The Dirac Theory — Pro and Con: The basic equation of this theory is 
e2 


<a (i — wit}. (1.1) 


o| bo 


FPext + Preact = MU 5 Freact = 


As is immediately evident from a consideration of the homogeneous case (Fex; = 0) 
this equation may contain physically unacceptable solutions. These are the so 
called ,,runaways“ in which the electron self-accelerates to infinite energies. The 
occurrence of these non-physical solutions has given rise to a long history of 
controversy both about the domain of validity of (1.1), see for example [86], and 
also regarding the possibility of seeking subsidiary ,,regularizing*‘ conditions 
that would automatically sort out the physical solutions. At the present time these 
controversies seem to have come to a stand-off: Every fresh objection or excep- 
tion to the application of (1.1) has been met with a severer regularization method 
or some other stricture on the solutions. The question is therefore not so much 
one of isolating a mathematically well defined or physically acceptable sub-set of 
(1.1) but rather one of evaluating how rich the structure of this remaining theory 
actually is. In Section 4-B we introduce a general dynamical correspondence 
principle which permits a quantitative answer to this question. Given a particular 
Fx it is immediately possible to evaluate the deviation of the radiation reaction“ 
trajectory [i. e. (1.1) including Fyeact] from the trajectory with Freact = 0. The 
extent of the permissible deviation is then directly related to the richness of the 
structure of the theory. Various specialized criteria yield all previously known 
regularizing prescriptions. aha 

This dynamic comparison method is augmented by a general radiation 
condition in Section 5. This condition enforces the requirement that the radiation 
reaction work must (on the average) represent a transfer of energy from the 
particle to the radiation field. Analytically this condition can be stated in terms 
of a Fowler-Emden equation. A study of the characteristics of the solutions 
substantially verifies an old conjecture of Abraham that a proper energy transfer 
occurs only when the particle trajectories have a quasi-oscillatory component. 


(ii) Runaway Solutions in Non-Local Electrodynamics: A number of covariant 
linear generalizations of the Maxwell-Lorentz theory have been proposed in 
which invariant form factors are introduced to eliminate the self-energy diver- 
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gences. (These methods are the precursors of the renormalization program in 
quantum electrodynamics [Bopp, STUCKELBERG, McManvs, FryNMAN].) It was 
first pointed out by LeHMANN that all of these form factor theories are essentially 
equivalent to a regularization (in the sense of Pauit and VILLaRs) of ordinary 
electrodynamics with a continuum of negative energy meson fields. As a conse- 
quence the electrons in these theories all acquire a negative non-electromagnetic 
mass. A general result of Hercguorz and WitpERMUTH then shows that these 
form factor theories necessarily must contain runaway solutions. Similar pro- 
perties have also been verified for a truncated version of quantum electrody- 
namics by VAN Kampen and Norron and Watson. 
The original extended electron theories of ABRAHAM, LORENTZ, and SOMMERFELD 
on the other hand do not contain any runaway solutions. The mathematical 
reasons underlying this distinction are discussed in some detail in Section 6-F. 
It is emphasized that the occurrence of runaway difficulties in the form factor 
theories reopens the entire question of a proper relativistic generalization of the 
extended electron concept. 


(iii) ‘Collective Modes and Excited Charge States: Electrodynamics in principle 
contains a “strong coupling” limit corresponding to the possibility of processes 


dominated by the closed loop cycle CQL (see the block diagram above). These © 


cycles represent transient collective modes in which radiative losses are greatly 
damped relative to their usual rates. Point particle theories of the kind represented 
by (1.1) do not contain a description of these strong coupling phenomena: The 
physical solutions in these theories are confined to two types of high frequency — 
behavior — (1) either nothing at all happens as the frequencies are increased past 
the critical frequency ~mc?/e?; or else (2) the solutions cease to exist altogether 
beyond this threshold. In order to obtain a non-trivial description of these col- 
lective modes it is necessary to introduce some kind of non-local electrodynamics. 
‘(This situation has some similarities to the “superconductor solution” problem 
in quantized field theories [NamBu, GoLDSTONE].) 

The basic equations of the non-local theories [e. g. (6.3) or (6.12b) of Section 6] 
show that these quasi-radiationless states can exist only at certain characteristic 
frequencies and with definite times of decay. This electromagnetic ‘‘isobar’’ 
spectrum is completely determined by the details of the particle form factors. 
‘Both classical and quantum mechanical estimates indicate that the ratio of the 
electromagnetic isobar energy to the ground state energy may be surprisingly 
large — of the order of 102 or greater. 

Following a suggestion of Boum and WEINSTEIN we discuss the implications of 
this theory for the elementary particle problem. In particular we present an electron 
isobar formula which predicts a mass value close to that of the u-meson. 


(iv) The Phenomenological Theory of Radiation Reaction (GrxzBuRG and Exp- 
MAN): The dynamical effects of radiation reaction may become significant in 
connection with certain coherent radiative loss mechanisms in bulk matter: 
Cases of practical interest include the Cerenkov effect — in cavities as well as in 
matter; transition radiation, Cerenkov self- excitation, and the complex Doppler 
effect. The object of the Ginzburg-Eidman theory is to provide a unified des- 
cription of the dynamical aspects in all these cases. In a sense this theory may be 
considered to be a macroscopic ‘“‘boundary condition” on all the microscopic 
theories discussed previously. No real link between any of these theories has 
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however as yet been established. A point which may be of interest for future 
work is that the mathematical structure of the Ginzburg-Hidman theory most 
closely resembles that of the original Sommerfeld extended electron theory. 


2. Radiation Reaction Theories Based on Conservation Laws 


A. STEWART and LarRMoR 


The beginnings of radiation reaction go back at least as far as 1871 when BALFouR 
Stewart [Z], on the basis of qualitative thermodynamic arguments, came to 
the conclusion that a moving body must be subject to a retardation owing to 
its own radiation. Sir JosepH Larmor took up the subject in an extended series 
of memoirs running from 1895 to 1927. Originally basing himself on a detailed 
aether theory and later relying on more general conservation arguments, he 
concluded that there should be a dissipative force accompanying radiation given by 


te 


Pent =e Pe Fol v, 


(2.1) 


where F is the rate of radiation and v is velocity.) In its final form the reasoning 
was that since in general total force is given by 


d dv dm 
Fic = di (mv) =m aa + © di’ 


and since in the case of a radiating body it is plausible that dm/dt =— R/c?, this 
“extra’’ force component ought to be identified with a radiation reaction. LARMOR’S 
formula will provide an interesting contrast in a subsequent discussion (Section 


3-C). 
B. The Radiating Oscillator 


The modern work was initiated by PLANcK [3], ABRAHAM [4], and Lorentz [5] 
with an analysis of the energy balance for a radiating oscillator. From the Lién- 
ard-Wiechert potentials it follows quite generally that the energy radiated in 


time 7’ is 
Vy 


t+ 
Dee 
cites { #2 dt. (2.2) 
Cc 
t 


The oscillator energy H is given by 
; Basho 4 Kha. 


The radiated energy loss is then presumed. to be compensated by a decrease in H. 
Balancing the energy flow yields the requirement 
aoory 


dB 2@.. f 
ae cage | aap (2.3) 


t 


1) See his Collected Works [2]; especially vol. I, pp. 414—597, and vol. IT, pp. 420—449. 
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One cannot, of course, from this alone draw any unique conclusions about the 
details of the oscillator dynamics. Nevertheless it is possible to get an equation of 
motion simply by setting the integrand of (2.3) equal to zero. This yields 


Uicgmd 5 Kite es (2.4) 


which is a non-linear equation. At the cost of some approximation it is also pos- 
sible to obtain a linear equation. Suppose that (2.3) were rewritten in the form 


Sipe Dee PAE 

Cores Ce ae 
Le (2 eet ii) — Ga} t= 0, 
t 


Then imposing the restriction 


5 me? 
(K/m)'? < a 1023 sec“, 
it follows that 
: DR A ee 


Finally, neglecting this last term in the integrand, we obtain 


2 
me =" — (oF a SE Vals A 
mx = —Ka+ 3 poe. (2.5) 


and this is a linear equation. (The consequences of the limit interchange are 
discussed in Section 3-B.) 

We should like to take some care in exhibiting still another approach to this 
simple situation: Supposing one insisted from the outset that there were some 
force Freact Which when added to F = mv would reproduce the radiation losses. 
Instead of (2.3) we should then write 


ea? t+7 t+7 
ee t.. 
= q)®e- — { Pras dt, (2.6) 


t t 


Integrating by parts on the left and choosing 7’ such that % vanishes at each end 
of the interval, there results 


Freact = —>- oe x (2.7) 
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to obtain a solution. Note that this is consistent with the previous result (2.4 
_ but with the important difference that the frequency restriction (K/m)' < 10% 
has fallen away. 


Since the details of the nature of the oscillator do not seem to be involved in this 


argument, it is now inviting to assert that (2.7) is in fact a generally valid formula. 
This assertion may be formally stated in terms of two axioms: 


(I) The dynamical effects of radiation reaction may be taken into account through 
an augmented form of Newton’s second law. 


(IL) One universal form of Freact is sufficient for all cases. 


The entire radiation reaction theory may then be summarized in terms of the single 
equation 


24 62... s 
Fext + Tee eee (2.8) 


C. Relativistic Generalizations 


Still another reason for selecting (2.8) as the basic equation of a radiation reaction 
theory is its natural place in a relativistic setting. It is well known, for example, 
that when »/c is no longer a negligible quantity (2.2) must be replaced by 


2 e* is Bigs es 

See BE aes * 

3 aay fe +|2 vil hat (2.9) 
t 


where y = (1 — v?/c?)-‘/. If one now retraces the steps leading from (2.6) to 
(2.7), it may be shown that the complete (vector) form of the result is (ABRAHAM[4]) 


Fna= = G7 {6 +(L) @- He +3(Z)@-H o+3(Z)e@-arel 
(2.10) 


and this is exactly what would have been obtained by a straightforward Lorentz 
transformation of (2.7) (von Lav [6]). With the help of (2.10) it is then a simple 
matter to transcribe (2.8) into its fully covariant form. We obtain 


2 e (d?u% du;\? du, 

where F,, and wu, are the Minkowski force and four velocity respectively. It should 
be noted that this involves more than the simple replacement x —> Bulds*; a 
point which in more formal treatments (e. g. [7]) requires special discussion. 
Tn the case of one dimensional motion, relativistic modifications do not bring in 
any essentially new features. On very general grounds [8] it is to be expected that 
the change of variable 


v = csinh (w/c) (2.12) 
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formally reduces these equations to the non-relativistic case. Indeed substituting 
(2.12) into (2.11) we find that : 


es Salo : 
Fext + "a cb ah te ee 


i.e. formal coincidence with the non-relativistic version except that now the 
proper time (s) has assumed the role of the independent variable.”) 


D. Dirac’s Subtraction Formalism 


The preceding considerations have been largely phenomenological. No direct 
attempt has been made to link the resulting radiation reaction equations to the 
Maxwell-Lorentz system; and it is in fact not even clear whether they are 
contained within this system or whether they involve additional assumptions. 
In 1938, Dirac [11, 12], stimulated by the then prevailing divergence difficulties 
in quantum electrodynamics, attempted the construction of a relativistic radiation 
reaction theory from first principles. [The efforts to formulate a satisfactory 
classical model theory have continued since that time (e. g. [13], [14]).] There are 
two main points involved in Drrac’s approach: 


(i) A mixing of retarded and advanced potentials to eliminate the diver- 
gences of the point electron theory. 

(ii) A more careful derivation of the equations of motion from the conser- 
vation laws by means of Gauss’ theorem. 


We shall discuss these separately. 
(i) Supposing one wanted the most direct approach to a radiation reaction 
equation. Within the Maxwell-Lorentz system the starting point then 
must be 


eF,; U; = m= (2.14) 


since it is the only equation of motion available. The electromagnetic field tensor 
may be thought of as being split into two components 


Fy, = Fett Peat, (2.15) 


i. e., an external field created by distant sources and a self-field due to the charge 
in question. For F<! we may in turn make the ansatz 


Feet =e Fee Fred, (2.15a 
r] 


where F';% represents the “convected’’ Coulomb field and Fi" denotes that part 
of the field which asymptotically appears as radiation. At the position of the 
electron F'%' of course diverges and it would be meaningless to substitute (2.15) 
and (2.15a) back into (2.14). However we certainly know that the Coulomb part 

i; has a divergence; in this case there is the hope that this is the only place 
where the infinities of F}!' appear and that F%* at the position of the charge 


My 6, reductions goes back to WussEt [9]. Extensive applications may be found in Puass 
0). 


SL ee, 
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| remains finite. Since it is presumably only this additional radiation component 


which affects the motion, we may simply omit Ff, and tentatively assume 


eeu; + eFPRdy, = m aes (2.15b) 
8 
as the final equation. In this way the radiation reaction effects emerge as a direct 
and plausible consequence of the Lorentz equation (2.14). 
A number of heuristic arguments are now available for evaluating F734, For 
instance, if Fj/' is assumed to be a real valued function of time with an appro- 


priate range of definition than it is always possible to decompose it into even and 
odd parts, viz. 


Pet = 4 (F(+) + F() +4 sth) — Pet}. (2.150) 


Fis (+) corresponds to a retarded solution of Maxwell’s equations and 
Fis! (—) corresponds to an advanced solution. The even part of F¥elt is then 
symmetric under time reversal and does not contribute to the radiation field 
(Pace [15], Norpstrom [/6]). The odd part however does change its sign and this 
is a sufficient condition for describing energy dissipation.’) Identifying this term 
with F724 and evaluating it at the position of the electron indeed gives a finite 
result; in fact the very same result which already stands in the previous relativistic 
formula (2.11). 

So far, so consistent. It should however be cautioned that the omission of Ff and 
the various field decompositions are heuristic devices which require careful justifica- 
tion (Haaq [18]) and are by no means the unique possibilities. [See for example 
ELIEZER [19], BERGMANN [20], and SCHONBERG [21] for variant theories.] A more 
systematic method of obtaining the equations of motion, without the Coulomb 
divergences, can be based on the conservation laws: 

(ii) Consider a cylindrical region 2 of space time pierced by a thin tube of (in- 
variant) radius r containing the world line of an electron. Assuming 2 to be 
source free, the divergence of the energy-momentum tensor vanishes in this 
region, and by Gauss’ theorem 


(ee 42 = § Ti do, = 0 (2.16) 


OX, 


' where o, denotes the surface of Q. Neglecting the contributions from remote 


space-like regions, (2.16) simply equates the energy-momentum flux across the 
tube surface with the net flow at the time-like ends of Q. The tube contribution 
may be expessed in terms of a single line integral extended along its length 
1 C21; 
xt oh eS oea - 241 
| fen, [re + Fe] —=—*+ O(n} ds, (2.16a) 
where the second (divergent) term represents the flux due to the Coulomb field 
F Es. If we impose the condition that the integrand be a perfect differential, i. e., 
1 eu . 
eu leer a FA Sees aE + O(r) = B, (2.16b) 
3) On this point see ScHwiNGER [17]. The dissipative nature of the radiation reaction does not 
however introduce any intrinsic irreversibility into electrodynamics; see PLANCK [3]. 
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for some four vector B,, then (2.16a) depends only on the flux at the ends of Q. 
However just as in the previous discussions (Section 2—B) this is not sufficient to 
determine an equation of motion since the conservation laws do not determine 
a unique choice of B, (BHaBHa [22]). A plausible selection which has the virtue of 
simplicity is 


Res (m ai aes “) Up. (2.16) 


This makes possible the transition to a point electron (lim r —> 0) on both sides 
of (2.16b) and yields an equation of motion which once again coincides with (2.11). 
It is apparent that this procedure makes the F£¢ subtraction less arbitrary and 
also shows clearly that it has the character of a mass renomalization. An important 
feature of this approach is that since no characteristic thresholds appear in the 
derivation, there is hope that (2.11) can in a sense be considered an exact equation.*) 


3. Discussion of the Solutions 


A. Physical and Non- Physical Solutions 


In the preceding section we have listed a number of arguments which converge 


on the expression 
2 


Pa eh as x 
Fext + 55 ey = MX (3.1) 


and its relativistic analogue as being the exact classical equations of motion 
including radiation reaction. In the present section we shall not present any further 
justifications of (3.1), but rather accept it as the final statement of a theory and 
turn to the consideration of various predictions and consequences which follow 
from it. 

We first remark that due to the general feebleness of the effects of radiation 
reaction (cf. the remarks in the Introduction) it is to be expected that solutions of 
(3.1) should be very close to those of 


Fext = Mx, (3.1a) 


i.e. the usual equation of motion omitting Fyeac¢. However, the addition of 
an x term to (3.1a) not only perturbs its solutions, it also adds entirely new 
solutions which may or may not have a physical significance. In order to see this 
more clearly, let us assume for the moment that the Laplace transforms of (3.1) 
and (3.1a) exist. Introducing the notation 


for the transforms; and setting 2e?/(3mc’) = + we have for the Laplace trans- 
form of (3.1) 


1 oe 
— Fox + pic! = pi’. (3.2) 


4) This was in fact Dirac’s original hope. For his second thoughts on the matter see [23]. 
Generalizations of the ansatz (2.16c) have been discussed by ELIEZER [24]. See also Section 5. 
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ima ext i t 


(278) ; (3.2 a) 


Solving for p, we find 


p= 


and this clearly shows the existence of the extra solution. In the limit + —> 0, the 
solution with the minus sign goes over into the equivalent of (3.1 a) and therefore 
may be considered as the “physical” solution; the plus sign yields a solution that 
diverges as t—>0 and represents a parasitic or ‘runaway’ solution. This con- 
venient identification fails to extend to the general case however. If for instance 
F ext is specified as a function of x and not of ft, then (3.2) must be understood as 


an integral equation for x(t). Fé; may then acquire an explicit + dependence. 
In this event 


lim {- sae G! Fe (3.2b) 
m 


™>0 


is no longer necessarily zero and neither solution of (3.2a) is constrained to reduce 
to a solution of (3.1a). These and related features may best be discussed with the 
help of a number of specific examples.*) 


B. The Damped Harmonie Oscillator 


This problem has been considered in some detail by Lorentz [5], Puancx [3], 


Puass [10], and Lorncer [26]. If we set Fox = —mw2x and continue with the 
abbreviation +t = 2e?/3mc® then (3.1) becomes 
tT —%—olx =O. (3.3) 


The general solution may be written as 
3 
Sead Oech (3.3a) 
=1 


where the /,’s. are the roots of the characteristic equation 
th —-P—owj =0. (3.3b) 


This cubic always has one positive real root and two conjugate complex roots 
with negative real parts. For the real root, we have 


Trae Aad) (3.30) - 


where it has been convenient to introduce the abbreviations 


5 if, 
al -les + VE ! ar aty| (3.34) 
0 Yo 


It is easily seen that in both of the limiting cases oy < ies (Sh orice and 
Wy >t (=> l, ~ [o,2/7]""); 1, represents a catastrophically rapid “runaway 


5) From the mathematical point of view, these difficulties stem from the non-uniformity of 
the 1+->0 limit. A rigorous discussion of this aspect of the Dirac equation may be 
found in [87]. Formally the situation is closely analogous to the small viscosity limit 
of the Navier-Stokes equation [25]. 
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solution. [This illustrates the remarks following (3.2a).] Imposing now the initial 
conditions (0) = 0 and £(0) =v, it is possible to suppress the runaway by 
choosing C, = 0. The complete “physically reasonable” solution of (3.3) then 
may be written in the form 


a =. e-*tsin wt (3.3) 
@ 


with 


3 0} 
Bee: ade (42 — B?) 


where A and B are as given in (3.3d) above. It should be noted that all the con- 
stants of integration have now been exhausted; the value of 2 (0) is automatically 
fixed at —2av, by the choice of the other initial conditions including also the eli- 
mination of the e4! solution. 

We shall return to this point later. 

The extremes of large and small w, produce some simplifications in (3.3d) and 
(3.3e); we find: 


Oy ts Oy tas 


(3.3f) 


A remarkable feature of these solutions, which can easily be read off the structure 
in the extreme cases, is that despite the presence of a natural threshold in this 
problem, i.e., t+, the nature of the solutions does not change even for the highest 
(@) >t) values of wy. This is in marked contrast to the behavior in other cases, 
e.g. the potential step (Section 3-D), where a characteristic threshold does make 
itself felt and plays the role of an upper limiting cut-off for the validity of the 
classical theory. However the structure of the cubic (3.3b) is sufficiently trivial so 
that neither “‘critical’’ frequencies, nor fields, or distances can affect the behavior 
of the oscillator. Moreover one can check directly that the details of the energy 
balance 
fo2) 


cm [wade = — Em (3.32) 
0 


wes 


are completely insensitive to the values assumed by «@,.°) 

It is interesting to contrast this with the behavior of the alternative ‘energy’ 
equation (2.4) which also admits of solutions of this type. In this case the character- 
istic equation corresponding to (3.3b) becomes 


tB + 7 + w= 0. 


6) This elegant treatment of the oscillator is due to A. CosTrKAs. 


/ 
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This cubic has a negative real root and two conjugate complex roots with positive 
real parts. In the w)<t limit 
; TOs 


lL, = —r-; ly ee F iw, 


so that the solutions represent either a tremendously overdamped decay or a 
slowy increasing oscillation. 


C. Time-Dependent External Forces 


If Fx is specified explicitly as a function of time, then the motion may be 
obtained from (3.1) formally by quadratures. The general result is 


@ = a + uot + at? fe a ( as =} Da 
T 


t t ts 
ellt 
_ at ff et Hen (tz) dt, dt, di, (3.4) 
0.0 0 


where 2, vp, and a, denote the position, velocity, and acceleration at ¢ = 0 
respectively. The ‘runaway’ factors e4/* signalize the presence of the unphysical 
solutions discussed previously [see (3.2a)]. In analogy to the treatment of the 
oscillator we shall see that it is possible in general to separate out the dynamically 
acceptable trajectories by proper adjustment of the initial conditions. Consider 
first the simple case of a step function force 


O aioe ged 
reali ator as 10, 


(3.5) 


From the general formula (3.4) one can immediately read off the solutions 
(note x = 0), 


F 
pee. Ore ge = vt} Po {2 + Qtr 4+ 22} 4+ zetlt Ja, — 81 —agr(r —#) (85a) 
— 2m m 


t 
(<0; m= 0,t-+-¢,t* ie" — (: + =f (3.5b) 
Clearly (3.5a) has runaway components. From the structure of the solution it is 
however evident that these may be cancelled out through a proper choice of ap, 
that is a) = Fy/m. With this choice and the auxiliary condition «(— co) = 0, 
~ the dynamically reasonable portion of (3.5) may be rewritten as 


iF, : 
eee t>0 3.5¢ 
vot + 5B; > (3.5¢) 
Ho (ett — 1}; #<0. (3.54) 
m 


From a physical point of view this implies that we have managed to get rid of the 
divergent solutions at the cost of introducing a slight acausality in the electron’s 


356 | ‘THomas ERBER 


motion. Since this electron ‘pre-acceleration’ involves only distances comparable 
to the electron radius (or times of the order r) it is actually more of a conceptual 
departure than an experimental difficulty. In fact this ties in very nicely with the — 
elimination of the “extra” boundary condition (i. e., a») introduced by the x term 
in (3.1). The proper ‘phasing out’ of the diverging solutions provides a natural 
employment for this otherwise superfluous element. 

Although this procedure is dynamically acceptable, it still does not yield a com- 
pletely reasonable solution. The root of the trouble lies in (3.5¢): Not only has the 
adjustment of a) removed the runaways, it has in fact removed all the traces of | 
the radiation damping. The trajectory (3.5¢) is completely identical to that given 

by the usual form of Newton’s laws (3.1a). We have ~=0 for ¢>0, and 
consequently the work of radiation reaction vanishes identically. If we simul- 

taneously accept the fact that the accelerated charge radiates then it seems as 
if the energy balance (2.6) no longer describes the situation. 

This is a non-trivial problem, having some bearing on the equivalence principle’), 

which has been discussed extensively in the literature (FULTON and RourRuicu [28], 

Rouruicu [88] and further references cited therein). There are a number of 
authors (e. g. DrucKEY [29]) who have attempted to rectify this situation by an 
artifice resting on an integration by parts: i. e., 


ty + 
. ae dt. (3.6) 


Com 
t 


ty 
| Breast edt=xi 
ty ; 


1 


The interpretation attached to this equation is the following: The first term on 
the right hand side is supposed to represent an influx of energy from the field in 
the vicinity of the particle—the so called “acceleration” or ScHott [30] energy — 
which then reappears in the second term as radiation to the far field zone of the 
particle. Since the two terms balance, Fyeac, may be zero although there is still 
radiation. This interpretation is however clearly contrary to the essential spirit 
of the radiation reaction development to this point. We have consistently sought — 
to identify the origin of radiated energy in the work done by the particle on the 
field: In the ‘acceleration energy’ argument the accelerated particle becomes me- 
rely some kind of transducer which transforms near field energy into far field 
energy. It should be noted further more that even the most detailed studies of 
the structure of the electromagnetic field in the case of uniform acceleration 
(SomMERFELD [31], LinpEMANN [32]) have failed to turn up any trace of this 
kind of behavior.’) 

This interpretation becomes even less tenable if we consider a slight modification 


of (3.5), viz., 
(torent 0) 
Post = 
Bot sfor et 10 (3.7) 


In this case the solutions still show signs of radiation reaction, even after the 
excision of the runaways, but the modifications are in the wrong direction. One 


may easily verify that iF Freact * dt has now become a positive quantity so that 


”) B.S. pe Wirr and R. W. Breunme [27]; F. RoHRLICcH (unpublished) ; S. CoLemaN, private 
communication. 


*) This point has also been emphasized by Rongxice, [88] and private communication. 
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in effect there is work done on the particle by the field; the kinetic energy is 
increased at the expense of the field in excess of the work supplied by (3.7). 
There is of course another way of meeting this problem, and that is to simply 
agree that Preacs = 0 implies no radiation at all. (See for example [10], p. 27.) 
In the context of the uniform acceleration problem this may appear to be an 
absurd conclusion but it certainly is a consistent consequence of the radiation 
reaction axioms listed in Section 2-B. It is in fact precisely the conclusion that has 
prevailed for years with respect to the hyperbolic motion of special relativity. 
(von Lave [33], Born [34], Pautt [35]) If one calculates the relativistic form of 
the radiation reaction (2.10) for the hyperbolic motion one eventually finds that 
Freact = 0. From this it has often been concluded that the hyperbolic motion 
does not give rise to radiation. From our point of view it is evident that these two 
problems, the relativistic and the non-relativistic constant force, are simply 
related by the hyperbolic transformation (2.12). It then becomes clear that such 
“proofs” of the radiationless character of the hyperbolic motion are about as sound 
_ as the corresponding claims for the non-relativistic case. FULTON and RowRuicH [28] 
have recently given an extensive summary of this situation and by a direct recal- 
culation of the fields, using a covariant generalization of the Poynting vector, have 
finally confirmed the existence of radiation. The “no radiation’ hypothesis is 
therefore untenable and the physical acceptability of (3.5¢ and d) is in doubt. 
By contrast, the Larmor formulation (2.1) seems to do very well in providing 
reasonable trajectories for (3.5): If for example the radiation rate R is assumed 
constant then the particle will simply follow a velocity damped trajectory. In the 
more realistic case that R is identified with x, i. e. (2.2), then (2.1) and (1.2) lead 
to the non-linear equation 

1 


(ca tiliee we 
ja ee (3.8) 


The variables may be separated with the result 


dv C tly |% 
a 4. — : 
dx 2rv? {1 aa 1 ah me 4 
Integrating, using the initial condition x = 0=> v = vq; and introducing 


« = 4 Fyt/mc*, one finds 


pl {F F aig Baw — 2)(1 + ao) + (terms in wah Gers: 


In order to identify the physical solutions we consider again the limit +— 0. 
- In this case (3.8a) reduces to 


Dye v2? — v3, 


and evidently the lower (++) sign is the appropriate choice. In the ‘strong’ radia- 
tion limit «av > 1 (3.8a) becomes 


ah 
a ld 
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which clearly exhibits the radiation damping in contrast to the usual v ails 
behavior. From (3.8) it also follows that ultimately R ~ a~‘/s, i.e., the radiation 
losses slowly decrease. 


D. The Potential Step 


Our examples so far have illustrated various dynamical difficulties associated 
with the appearance of divergent runaway solutions as well as energy difficulties 
associated with the behavior of the radiation reaction work. There are still other 
difficulties latent in the % feature of (3.1) as elegantly exhibited in an example 
originally due to Bopp [13]. 49 
Consider a potential step of height V, situated at «= 0. If we integrate the 

equation of motion 


2a 7 ees (3.9) 


across the potential step from left to right and assume the continuity of x and 2, 
then we find the condition that 


Vo 


MT Vo 


a (0+) — a (0-) = 


where a,(0+) represents the accelerations immediately before and after x = 0, 
and v, = «(0). Let «(0) = 0 and designate the initial velocity of the particle 
by v;[= x2(—oo)]. With these conditions (3.4) may be rewritten as 


a= vt + t(vp — v,) (e*® — 1); t<0 (3.9a) 


v= (u— Te) t+r(y—m + 22] (e/= —1); §>0. (3.9b) 


If the particle is to penetrate the potential step and continue towards the right, 
the runaway component of (3.9b) must be eliminated. From (3.9b) it becomes 
clear however that this requirement will not specify a unique trajectory. Setting 


Vv 
Dordt AO 
0 


in order to suppress the exponential, we have 


i 4V,|2 
U9 = {! = 1 = | I, (3.10) 


a 


and therefore two values of vy are compatible with all the restrictions so far. 
In particular, if 4V,/mvj <1, it follows from (3.10) and (3.9b) that 


Se (3.11 a) 


(3.11) 


i 
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are both possible final velocities of the particle. In terms of the dynamical criteria 
advanced so far it is impossible to distinguish one of these as being more reason- 
able than the other. ei 

If we now draw on energy considerations (3.11 a) appears to be a somewhat more 
plausible choice: One can easily verify that the kinetic energy behind the potential 
step is diminished essentially by V, with a small additional amount (~ V,/mvi) 
being expended in radiation. (3.11b) on the other hand represents a motion in 
which practically all the energy is radiated away at the potential step.®) We 
note further that (3.10) imposes $ mv,?2 > 2 V, asa minimum energy requirement. 


In order to eliminate any artificialities that might be due to the abrupt rise of 
the potential, Haag [78] has extended this investigation to the case of a uniformly 
increasing potential. He finds the interesting circumstance, reminiscent of the 
Klein paradox, that as long as the potential rise is a gradual one over a distance 
of the order of cr, the duplicity of solutions may be avoided. Here again one en- 
counters the elusive feature of a critical distance, or ‘‘electron radius’’, in a for- 
mally point electron theory. WEssEt [9], StickeLBeEre [37] and Feynman [38] 
have in fact used this multiple solution property of (3.9) to discuss the possibility 
of a classical description of pair production. 

Another curious feature displayed by this example is that in the case of solutions 
of the type (3.11b) there is an invariant upper limit to the amount of kinetic 
energy that may be transmitted across a potential barrier. This follows readily 
from (3.10) and (3.11b) since together they imply the inequality a(t > 0) 
< [m/4V,]:. Pomprancuuk [39] has elaborated this argument into one of the 
few non-trivial applications of the radiation reaction theory. He showed that if 
the potential barrier is identified with the earth’s magnetic field, then in the ultra- 
relativistic case any electronic “‘primary”’ cosmic ray component would be limited 
to energies of less than 101” ev at the eath’s surface —independent of the magnitude 
of the initial energy. This point is also connected with the interesting possibility 
of a vacuum Cerenkov effect. 


E. Singular Cases 


The particulary interesting possibility that no solution of the radiation reaction 
equation resembles the solutions of the corresponding equation omitting Preact 
was first discussed by Etrezer [24]. From our previous remarks (see 3.2b e¢ seq) 
- it follows that this is equivalent to constructing an F’x4 (w) which when re-expressed. 
as a function of time acquires a singular t-dependence. EL1EzER showed that the 
one dimensional Coulomb force satisfies this condition and in particular that 
radiation reaction will cause two unlike changes to repel sufficiently so that no 
collision is ever possible. This highly paradoxical result has occasionally been 
dismissed on the grounds that the Coulomb potential is too singular to admit 
solutions near its origin. It has also been argued that quantum mechanical effects 
can be trusted to intervene at sufficiently small distances to perturb the classical 
idealization in the right direction. It is difficult to agree with these points of 
view for at least two reasons: iicdy 

(1) Despite the singular character of the Coulomb potential it-is only necessary 
to change it from an’ attractive to a repulsive potential in order to introduce 
solutions representing collisions (ZrN [40]). Solutions approaching the origin are 


\ 


2) ‘A detailed treatment of the relativistic 6-function force has been given by: AsHaumR [36]. 


26 Zeitschrift ,,Fortschritte der Physik’ 
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possible but just under those circumstances directly contrary to what one would 
ordinarily expect. 
(2) The use of a quantum mechanical argument at this stage would be an inversion 
of the logic of the development. The difficulties underlying the classical theories 
are to a large measure shared by the quantized versions. On a later occasion 
(Section 6-F) we shall in fact use some of the insight gained. from our examination 
of the classical theories to comment on various “runaway” features which also 
appear in quantum electrodynamics. (See especially footnote 21.) 

A number of other singular cases involving both bounded and unbounded poten- 
tials have been discussed by WILDERMUTH [47]. 


4, Physical Solutions: The General Case 


A. Regularization and Perturbation Methods 


From the preceding it has become clear that we cannot accept all solutions of the 
radiation reaction equations as physically meaningful; some further procedure 
must be adjoined to the integrations in order to generate the physically reasonable 
solutions. In this connection, a number of schemes have been proposed by Haaa, 
WILpERMUTH, Crnt, EtrezeR and others. We may group these efforts roughly 
as follows: (1) Regularization Procedures; and (2) Perturbation Methods. 

(1) Regularization Procedures: We have already seen by means of a number of 
examples that parasitic runaways are persistent features of (2.8) and (2.11). 
We have also seen that a proper adjustment of the initial acceleration can phase 
out the divergent components. Haag [15] and later Pass [10] have described a 
method which produces this phasing automatically. 

Consider (2.8) integrated once with respect to the time in the following way 


t 
1 
z — ellt ao = aa | oo Bnet . (4.1) 
0 


If in the limit ¢ > oo, %(t) is to remain finite it clearly is sufficient to choose 


1 
t= aa oth Fede. (4.1a) 
0 
Entering this back into (4.1) we may rewrite it in the compact form. 
zt = ae e-(—h)/t F t dt 
= — ext (h) ty, (4.2) 


t 


which displays the fact that the radiation reaction makes the particle acceleration 
proportional to the (incomplete) Laplace transform of the force — probing ahead 
an interval ~7z in time. This felicitous interpretation has encouraged some 
authors to advance (4.2) as the fundamental equation of the theory. *) As yet, 


a5) The relativistic generalization of (4.2) is extensively discussed in [88]; some related ques- 
tions regarding the existence of solutions are examined in [87]. 
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however, no physical argument has been given which would lead to such a result 
directly. It may be noted that the transition to (4.2) has an equivalent interpretation 
as a standard transformation in the theory of integral equations, i. e., the incor-. 
poration of boundary conditions by means of quadrature. It should not be over- 
looked that such a regularization procedure copes with only a part of the diffi- 
culties encountered in constructing physically meaningful solutions. Similar 
remarks also apply to the Fourier transform method suggested by WILDERMUTH [41] 
and to the asymptotic regularization of Exrezer and MativacaANaM [42]. See 
also Gora [43]. 

(2) Perturbation Methods: This is in principle a less satisfactory approach than 
the regularization procedure since from the outset all other considerations are 
subordinated to the aim of manufacturing reasonable solutions for (3.1) through 
perturbations of (3.1a). In particular this carries the penalty that occasionally 
one will construct perturbation solutions which actually have no counterparts 
among the solutions of the exact problem. [See also the remarks following (3.2a).] 
The approximation schemes in such cases may in fact conceal important features 
of the underlying theory. 

In the construction of these solutions one has the option of carrying out the 
t-perturbations either in the differential equation itself [e. g. (2.8)] or in its 
solutions. STEINWEDEL [44] and Buasuna [45] have noted that the unwanted 
components are usually characterized by terms in e/* which are non-analytic in 
t. They have therefore suggested following the second method: The solutions are 
to be obtained exactly from (2.8) and then examined for analyticity about 7 = 0; 
all singular terms are then to be discarded. A slight variant of this has also been 
proposed by Crnt [46]. ARLEY [85] has however criticized these recipes as being 
excessively stringent since too many otherwise well behaved functions would be 
excluded. Our previously acceptable “‘physical’ solutions (3.5¢ and d) provide an 
example of this since they clearly could not survive the analyticity test: 

The perturbation solution here would actually converge to a discontinuous 
trajectory ! 

The other alternative, i. e. modifying directly the differential equation, is actually 
the simplest and most foolproof perturbation technique. The radiation reaction 
term mt is simply evaluated from the unperturbed equation Fex¢ = ma and 
the exact equation (2.8) replaced by 


De eb et I (4.3) 


(or its relativistic analogue). This automatically yields the perturbed trajectories 
and voids all infinities ab initio. Indeed, judging by some of the current literature, 
for example [47], it is in fact (4.3) and not the prickly (2.8) which is the basic 
equation of the theory. 

Field theoretical perturbation schemes have also been discussed by STEINWEDEL 
[14] and Pricoarne and Lear [53]. 


B. A Dynamical Correspondence Principle 


We shall now show how most of the points of the preceding discussion can be 
summarized in terms of a precise mathematical argument. Consider again 


m x, == Disa (4.4a) 
26* 
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and the radiation reaction equation 
— tmx, + Mx, = Fest, (4.4b) 


which from this point on we shall regard as its ‘““degenerescent”’ partner.?°) 
What we should really like to study is the extent of the correspondence between 
the solution-sets of these two equations. 

Let us introduce a deviation function h(t) defined by 


—— 


h(t) = %_(t) — x(t). (4.5) — 


It is also convenient to suppose that the initial conditions of (4.4a and b) have 
been matched at t = 0 by requiring 


a,(0) =a,(0) and (0) = %,(0). (4.5a) 


Introducing now the formal solutions of (4.4a and b) into the right hand side of 
(4.5), we find 
t 


1 t 1 
h(t) = — on text + a se = ( ae = se { ef—h)it FO, (t,) di,, (4.6) 
0 


where F®, denotes the (indefinite) repeated integral of F.x¢, and a, stands for 
x,(0). The information contained in this equation is actually better displayed 
by inverting it to give F’, in terms of h. From this point of view (4.6) is really an 
integral equation. Its solution may be written as 


t 


Oo see mf hwa —mh + may 


0 


2 fi 
> (4.7) 
and this we shall regard as the basic expression of the dynamic correspondence be- 
tween (4.4a and b). 

The principle advantage of reformulating the problem in this way is that it 
preserves the freedom as to what may be considered a reasonable deviation 
between the solutions of (4.4a) and (4.4b). In particular, with a given bound on 
h(t), it is possible to make a value judgment as to just how rich the structure of 
the theory actually is. To illustrate this point consider the extreme case 


\h(t)| << K forallt. (4.8) 


From (4.7) it follows immediately that 

Fa Ola Ke: (4.8a) 
At most therefore FQ, ~é. If we now assume in addition that ¢- Fox; is contin- 
uous, increasing, and positive for ¢ > 0, then by a standard TauBERTAN theorem 


*) The terminology is adapted from the degeneration theory of differential equations 


(Mryorsxy [48]). 
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‘(e. g. [49]) we have the strong conclusion that 


PENG (4.8b) 


The uniform boundedness of h is obviously a very severe restriction.) 
Equation (4.7) also includes most of the preceding regularization schemes. If 
for instance Fx; = Fy, then (4.7) may be put into the form 


1 
h' — ro h = (ma, — Fy) t (4.9) 
which integrates to give 


h =t (may) — F,) {et*— 1— + (4.9a) 


Tt 


and clearly unless a = Fy/m, that is, h =0, there will be a runaway deviation. 
The same result could also have been obtained directly from (4.6). 
A more general situation is 


Pet (== 0 for O< t< tux, and. F'oxn(t) = 0. otherwise: 


i.e. an external force acting during a finite time interval. In this case we have 


t 
1 1 
h’ —-—h=a yt — — Be. , 
: ay 7 [ Fat a (4.9b) 
0 
The complete integration yields 


tmax tmax 


1 
h(t) = = | Foxld dt —a,1? + rel" Ja, — a fem Fext(t) dth.  (4.9¢) 
0 0 


In order to keep the deviation within reasonable bounds, clearly a) must be 
adjusted to 


tmax 


1 . 
Gg af Fext(t) dt, (4.9d) 
mT 
0 

and in the limit tmax — co this coincides exactly with the previous recipe, i.e. 
Haaq’s regularization. é 

This can be contrasted with the corresponding comparison for the perturbation 
equation (4.3): In this case it is easy to see that the deviation function is given by 


Aa = [Peal dt. 


Clearly, for reasonable Fx, runaway deviations such as those occurring in (4.7) 
are avoided. In particular, as t + 0, the deviation disappears, as is consistent 
with the perturbation approach. 


11) Analogous arguments can be carried through for weaker versions of (4.8). 


} 
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C. Instability of Solutions 


The solutions of the radiation reaction equation (2.8) have a further undesirable 
dynamical property and that is that none of them are stable! 

In the case of linear forcing functions this is equivalent to a well known result 
which goes back to the classical Routh-Hurwitz criteria. The general situation 
has recently been discussed by EzrrLo [50]. This result can be understood very 
simply in terms of our previous discussion: The key point is the delicacy of the 
adjustment needed to phase out the runaway solutions. If for example we consider 
a particular regularized solution corresponding to a given Fext, we can study the 
effect of a small additional perturbing force F@) which acts over t; Ct < fp. 
By a simple extension of the previous arguments it can readily be shown that the 
deviation from the regularized solution is of the order of 


ty 


efit 40 e-tlt F®) dt, 


for t > t,. In other words, small perturbations are impossible — all solutions are 
unstable. 


5. The Radiation Condition 


A. Introduction 


An acceptable radiation reaction theory must give a reasonable account of the 
transfer of energy between the particle and the radiation field. We have already 
seen that this is a decisive point in the case of uniformly accelerated motion: 
Although the radiation reaction equation (2.8) gives a reasonable description of 
the dynamics in this instance (once the regularization has been included), the 
work of the reaction forces fails to account for the radiation losses in a natural 
way. In gauging the richness of the theory it is therefore essential to take both 
dynamic and energetic criteria into account. In the present section we shall give 
a quantitative formulation of a radiation condition which will assure that the work 
of radiation reaction represents, on the average, a loss of energy by the particle to 
the radiation field. It will appear that this is a very severe limitation both on the 
admissible external forces as well as on the character of the particle trajectories. 
In particular we shall substantially confirm a conjecture originally due to ABRAHAM 
[4] that a consistent theory — based on (2.8) — is only possible in the case of quasi- 
oscillatory motions. 


B. Analytical Preliminaries 


It will facilitate the subsequent discussion if we introduce some standardized 
nomenclature. ; 

(a) Essentially increasing functions: A function f will be said to be an essentially 
increasing function on the interval Q if there exists a strictly increasing function 
g > 0 and a6 > 0 such that 


If) —9g@O| <0 (5.1) 
for all <Q. 
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(b) Quasi-periodic functions: Consider a function f, also defined on Q,. which has 


the property that there exists an ascending sequence of points ty, ty, . . .» try €Q 
such that (to) = f(t) =... = f(tn+,;). Assume further that there exists no other 
such sequence ,..., t, for which k > n + 1. The sequence {t;}1s then said to be a 


maximally repetitive sequence, Let 7; = ti, —t; be the repetition intervals 
corresponding to this maximal sequence. A quasi-periodic function is then defined 
by the following properties: 


(1) The partition {t,;} is a fine mesh of Q. 
(2) Max (t;) ~ Min (7;). (5.2) 
(3) f is uniformly bounded in Q. 


Conditions (1) and (2) essentially guarantee that the function repeats itself suf- 
ficiently often and with enough regularity to deserve being called quasi-periodic. 
(3) eliminates oscillatory function with arbitrarily large amplitudes of fluctuation. 
(c) Essentially increasing quasi-periodic functions: Combining (5.1) and (5.2) we 
now define an essentially increasing quasi-periodic function f in terms of the 
representation , 


f=ft+g. (5.3) 


Such a function is therefore the sum of a strictly increasing function and a quasi- 
periodic function. 

(It should be noted that the definitions (a) — (c) are not sufficiently precise to 
exclude a certain number of pathological cases [89]. Supplementary hypotheses will 
therefore have to be introduced several times in the subsequent discussion. In 
addition it will be convenient to assume that all functions which appear will be 
at least twice differentiable.) 


C. The Radiation Equation 


The total work of radiation reaction, in terms of the variables of (4.4b), is 
t 
W(t) = —mt ff tigty dt. (5.4) 
to 


On the average W should be an increasing function of the time: In general, of 
course, we cannot expect this to be a strict increase since the decomposition (3.6) 
suggests that there may be a fluctuating “acceleration energy” component 
representing a short term energy influx from the field. 
We therefore make the basic assumption that in any case W must be an essentially 
increasing function of the time. It is then an easy corollary of (5.1) that there must 
be some last point of time beyond which W (t) is certainly positive. Tt will turn 
out to be convenient to choose this point as the zero of time and to discuss (5.4) 
in terms of h(t) — the deviation function between the radiation reaction equation 
(4.4b) and the comparison problem (4.44) matched at this zero point of time. 
We begin by transposing (4.5) to give 
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Scaling x, (0) to 0, it follows that 


= 1 - ; 
pt Sako heer pay es (5.6) 
T T 
Introducing now the auxiliary variable 


ee peers (5.7) 
= 
(5.4) may be rewritten in the form 


t 
W(t) = —mrf iy dt (5.8) 
0 


where the fluctuating portions (if any) of W for t < 0 have been implicitly disre- 
garded.) In this form W is finally related directly to h(t) and it is possible to study 
the additional restrictions which arise as a consequence of W being an essentially 
increasing function. 
It is convenient to work with the differentiated form of (5.8). 

y + + a= 3 (5.9) 


This makes it clear that the radiation condition is actually equivalent to a modi- 
fied Fowler-Emden equation [51], [52]. The associated boundary conditions 
follow from (4.7) and (5.7): : 


1 : 
y(0) = 0; yO) = 45) 90) = G, = Mest). (5.10) 


In general, it is known that if W~O (t), m rational, then y ~ O (t*) for some k. 
More specific information will of course follow from our present restrictions on 
W and a number of other provisory assumptions. 

It will not be an essential restriction of generality to suppose that a, > 0. By the 
mean value theorem it then follows that there exists an interval Q, consisting of all 
t,0 <t< #t, for some t, > 0, where both y and y are positive. We first show that 
the simplest possibility, i.e. both h and W strictly increasing functions in Q,, 
cannot occur: It is convenient to introduce the equivalent conditions h(t) => 0; 


W(t) > 0, t« 29. Multiplying both sides of (5.9) by y and integrating, we have 


ite V 4 
“yay, (ht) eg Sy ee dt. (5.11) 


But by (5.7) and the present hypotheses, the left hand side of (5.11) is positive 
while the right hand side is negative — hence there is a contradiction and this 
simplest case is eliminated. 


1“) It will be convenient from here on to absorb the factor mt into W (é). 
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Suppose now one drops the condition on hf but still maintains W > 0. According 
_ to (5.9) y is then a concave function. This is a very severe restriction. It immediately 

implies for instance that no solutions of the form y = kt” exist; concavity restricts 
n to the interval 0 <n < 1; and this range in turn is eliminated by the finiteness 
_of 4 (0) (see 5.10). 
If y is concave then there is the possibility that y eventually becomes negative: 
From this it follows that there may be a ft, such that y itself becomes negative for 
all ¢ > t,. If however y(t.) is permitted to be equal to zero, this brings with it 
the consequence that W(t.) = 0 because (4.7) and the presumed finiteness of 
Fx imply that 7 cannot diverge. In this event a repetition of the argument leads 
to the conclusion that W itself acquires an oscillatory component — i. e. the very 
point we are trying to establish. 
Let us now go back for a moment and assume W > 0 and y> 0 in order to 
exhaust the alternative possibilities. It is shown in detail in Appendix 1 that under 
these circumstances 


ey = S yale 
yO) Sage }2 —exp| “—— [I (5.12) 
and also 
t 
| Jee 
t 
Wie) a2) 2-exp | = (5.13a) 
| aot 
Furthermore 
t 
IF Gi == (a2 ln 2) ¢; (5.13 b) 
0 


These inequalities constitute severe restrictions on the nature of W(t) and y(t). 
There is in fact a less explicit although still stronger result which may be stated 
in terms of the following ; 
Theorem: If y>0, then all (sufficiently differentiable) W functions compatible 
with (5.9) must satisfy W< t in the range t > 0.1) 

Proof: Integrate (5.9) by parts, 


t 
W Wy 
a — 9) => -/ z dt. (5.14) 
0 


The existence of the integral follows from (5.10). Since both y > 0 and W > 0, 
we have the inequality 


a — Y(t) = (5.15) 


ey Wat Saaue W/t = 0. See [49] for notation. 
— 00 
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Integrating once again, and introducing W®) — ‘if W dt, (5.15) becomes 
0 


O> YP —agty + WO. (5.16) 


In the second integration we require 


t 
(1) (1) 4 
us Oy and [= Gh, CS 
y (0) y 


which may be justified from the general assumptions about y and W. 
For a given W, (5.16) is an integral inequality for the permissible values of y. It 
may be rewritten in the form 


1] /2 (1) ]*/2 
x tlt = la }< y(t) < 2 gt f oe ae aH }. (5.17) 


From this it follows in turn that 
WY < 4 (apt). (5.18) 


The key point now is that one can also prove that WY ~ O(#) is impossible! 
From this statement the conclusions of the theorem will follow. 
Let 


W = at Hipiifor #> tee: 0 (5.19) 


and assume « > 0. Going over to the inverse functions, (5.9) is separable and 
integration yields 


y(t) {pg aunts SL ee 
ak barat a=) Latta) esse” cee O- (5.20) 
Therefore 


which clearly shows that y is bounded from above. On the other hand (5.17) and 
(5.19) imply that there is a lower limit on the growth of y which ultimately exceeds 
this value as ¢ — oo. We are forced to the conclusion that there exists a t; beyond 
which no solution to the problem exists ;!%) (5.19) must therefore be amended to 
read tj << ¢< ts, which finally implies that W cannot be ~ O(#). [If « and f are 
replaced by smooth bounded functions of ¢, the argument is essentially unchanged. ] 
Since a priori all pathologies have been excluded by (5.1) and (5.2) there are only 
the three possibilities remaining for W: Hither (1) W > #; (2) W ~¢; or (3) 
W <t. But (1) and (2) have been eliminated by the preceding; so this leaves 
only (3), and therefore completes the proof of the theorem! 

Summarizing all the work so far, one finally has the conclusion that if y is a 
non-decreasing function of ¢, then the radiation must proceed at a rate which is 
™) The equation may be completely integrated in terms of inverse error functions; t; then 
emerges from the bound of the domain of definition. 
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less than proportional to the time. In view of (2.2) this is clearly much too restric- 
tive a condition. Tt y is permitted to become negative, then either this reduces to 
the case previously discussed — through y eventually crossing the ¢ axis — or 


else fluctuations of y, reflected in jj, again force W to oscillate since jj and W are 
directly tied together by (5.9). The only other alternative is that both the deviation 
function as well as the radiation acquire oscillatory components. In this sense, with 
the exception of a number of highly circumscribed cases, we finally arrive back at 
ABRAHAM'S conjecture — which now may be regarded as confirmed — that from 


an energy transfer standpoint the theory yields satisfactory results only in the case 
of quasi-periodic motions. 


6. Non-Local and Extended Electron Theories 


A. Introduction 


A considerable portion of the theory of radiation reaction has been developed 
within the context of various non-local extensions of classical electrodynamics. 
Before proceeding with the detailed discussion of these features it will be useful to 
review some of the (contemporary) reasons for taking an interest in these develop- 
ments. 

(i) Non-Singular Form Factors and Runaways — A number of experimental 
programs are currently under way which are producing detailed information on 
the form factors of nucleons. If classical field theories, patterned after electro- 
dynamics, are to continue to be of any qualitative significance in discussing such 
particles then clearly these new structural degrees of freedom must somehow be 
accommodated. Inevitably this leads to some type of non-local generalization 
of electrodynamics. The associated radiation reaction theories then of course 
must also be superposed on some type of non-local structure. In the later deve- 
lopment it will appear that such additional constraints can actually be of help 
in making a proper choice among the many possible extensions of the classical 
theory. There are, for example, numerous “‘covariant form factor”’ generalizations 
of electrodynamics which can be shown to have non-physical radiation runaways 
despite the fact that the self-energy divergences have been successfully remo- 
ved. 

(ii). Phenomenological Theories of Radiation Reaction — Current progress in 
plasma physics and millimeter microwave generation has stimulated the 
development of an extensive phenomenological theory of radiation “reaction. 
(GinzBura and Erpman [54]) Although this theory has not as yet been put on a 
completely consistent foundation, it is already clear that its mathematical 
apparatus has only the vaguest resemblance to a point particle description. The 
spirit of the formulation is in fact very close to that of the classical extended 
electron self-force picture. an 

(iii) Electrodynamic Collective Modes — Non-local generalizations of elec- 
trodynamics are also interesting in that they offer certain formal possibilities 
which are completely lacking in the point charge theories. The most interesting 
of these features is the appearance — at certain critical frequencies — of self- 
sustaining charge oscillations. The charge oscillations have a natural interpre- 
tation as collective resonances of the field-matter system. In common with collec- 
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tive modes familiar from other areas of physics, these charge resonances are non-— 
analytic inthe coupling constant (e?) and cannot be reached in truncated versions 

of the theory. The existence of these collective modes offers the formal possibility — 
of extending electrodynamics into an “isobar” theory of excited charge states. 


B. Higher Derivative Theories 


The minimal intervention in the structure of the point charge theory which offers 
a way out of the contretemps of runaways and the non-existence of physical 
solutions, etc. is the introduction of higher derivatives. This avenue has been 
explored in detail by Exrezer. In his first attempt [79] Exrezer widened Dirac’s 
ansatz (2.16c) by including a new constant x, and some higher derivative terms as 


follows 
1 e 2 Puy | du; \? 
= SS G e al 
By=(m—5 \ ue txts at pu (Tf (6.1) 


If the original derivation is now repeated with this new form of B,, it is found that 
the resulting equations of motion contain fourth derivatives. The analysis of these 
equations becomes proportionately more difficult; however EL1ezER himself, and. 
later Zin [40] and Haag [18] showed that the runaway solutions and related 
difficulties still persist in this extended formulation. 

In a second attempt, ELrezEer [24] [55] introduced more drastic measures: The 
ansatz (2.16c¢) was generalized to include an infinite number of derivatives, viz. 


fie — 2 
B, =(m 9p) Me > Bon ue” — 


n=1 
on Neat t 
and E ue” ata = (— We use a) mS v+1) 25 eee 1)" au” us uf (6.2) 


where the B,, are arbitrary constants and u, is the four-velocity.In this form the 
theory is actually equivalent to an extended electron model.14) To see this, we 
introduce (6.2) into (2.16b): In the non-relativistic limit the result may be written 


m {3.4 >| re (6.2a) 
p= 


es ® represents the derivative operator and the c,’s are constants related to 
the B,,: 


In the special case c; = (t)//j! (6.2a) simplifies still further to 
ma(t+ t) = Foxt (6.2 b) 


by virtue of Taylor’s theorem. In this form the content of the theory is summa- 
rized in terms of a differential-difference equation; this in turn can be shown to 


coincide exactly with the corresponding expressions obtained from the extended 
electron theories. (See below.) 


™4) Note however that this approach still requires a mass renormalization. 


i] ) 
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C. Extended Electron Theories 
The original versions of electrodynamics were all of course based on the concept 
of an extended electron. As mentioned in the Introduction, the ultimate goal of 
these theories was the calculation of all of the inertial and dissipative properties 
_ of the electron from a model which pictured it as a sort of electromagnetic sponge: 
The various parts of this sponge were supposed to be capable of interacting with 
each other as well as with the external fields — the whole object responding to 
the net sum of the “outside” and “self” forces. For the fields, the basic equations 


were taken to be those of Maxwell. The dynamical properties were then assumed 
to follow from the equation 


Pest + Poor = Maneen © (6.3) 


which, because of the implicit a(¢) dependence of Fyeir, rather than being a 
differential equation is actually a functional equation.15) It is this functional 
equation property of (6.3) which endows the classical theory with its astonishing 
wealth of solutions! We shall see first of all that in the limit of small accelera- 
tions (6.3) reduces to (2.8), i.e. the ordinary theory of radiation damping. Se- 
condly we shall see that at “super-light”’ velocities (6.3) goes over naturally into 
the constant velocity radiation reaction of the Cerenkov effect. Finally it will 
appear that this equation also contains the radiationless charge self-oscillations 
discussed previously in Section 6-A (iii). 

Clearly this program requires a rigorous and generally valid expression for Fer. 
We begin by writing down the usual representation for the retarded potentials 


A(x’, t) =| ae" { aa” je — "8 flee) 0 (0" ! oe ]) (6.4) 


The Lorentz force is given by 


0A 
L(x’, t) = — grad {o _ 2 ze Al _ e{e x [x’ — x]) x curl A — a (6.5) 


where x locates the center of the electron and w is the angular velocity. The total 
self force is then 


Freir(t; 2) = f doe’ 2 (x, t) 0 (w',t) (6.6) 


where the x has been inserted on the left to serve as a reminder that Pep is a 


functional of the trajectory. wk) 
The sevenfold integration implied by (6.6) can be simplified in a number of ways. 
For example Fy.i¢ can be written as a Fourier transform 


Fseig = <a fdr fF lenox (ike) — x(t —t)]} X 
0 


x [kX (a(t — 1) X Kk) cos(ckt) + ika(t) X (a(t — t) Xk) sin(ckr)] (6.6a) 


15) The non-electromagnetic inertial component Mech is inserted as a purely amet Et: 
gical term. In a strict (classical) nucleon field theory the corresponding ansatz would be 


c i Yukawa field. 
Mmech & > F¥.14, where Fy; represents the influence of the 
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where ena 
on = (2n)-"lf g(a’, t) eth dae’ 


which reduces (6.6) to a fourfold integration. (See e. g. Boum and WEINSTEIN [56] 
or GrnzBurG and Erpman.[54].) We shall return to this point later in connection 
with the phenomenological theory of radiation reaction. _ ' ; 

The complete calculation of all the explicit integrals contained in this expression 
was carried out by SomMERFELD [3/] [57] for two special cases: The uniformly 


charged sphere and the infinitely thin spherical shell. For the uniformly charged 


sphere, in pure translation, the result may be written 


Ty 


fe oe (== g ier +9) — Siler +) ips 


T oT 


{i — x(t —T) a iler—T) aire (lex -7)) ath. (6.7) 
0 


The auxiliary function f is defined by 


z 
1) 364 12 8 
j= [ {ane 39+ 2 Fh ap= 5S P+ ent Ft OTe) 
0 


21o 


¢ denotes the path length traversed by the center of the electron in time rf, i. e. 
¢ = x(t) — a(t — 7). (6.7 b) 

Tt, and Tt, are the (possibly non-unique) roots of 
cr +7=2r, and |erv—( | =—2r, (6.7¢) 


respectively; ry is the charge radius. 

The functional nature of Fei¢ is now explicitly displayed on the right hand side 
of (6.7): This type of expression is actually an “‘integral-delay’’ representation 
because the retardations of the functions occuring in the integrand are themselves 
dependent on the unknown function. (See [58].) If this representation is carried 
back into (6.3), the resulting expression is evidently a functional equation in 
x(t). It should be emphasized that as far as the rigid spherical electron is concerned, 
this equation is an exact consequence of electrodynamics — presupposing no 
restrictions on either the velocity or the acceleration. 

It was of course one of the great triumphs of this classical theory that one could 
Set Mmech = 9, i. e. replace (6:3) by 


Pext erase (6.8) 
and so achieve a purely electromagnetic explanation of mass. (Pace the stresses 


of Poincaré.) These remarks depend on the use of the so-called ‘“‘quasi-stationary”’ 
approximation 


dv v 
pg aN eel: (6.9) 


The Classical Theories of Radiation Reaction 373 


which for ¢ — v reduces to the usual assumption that the accelerations are always 
small in time intervals of the order of the electron transit time (Za) systematic 
use of this inequality, i. e. expansions of the type ; 


a(t —t) ~ a(t) — a(t)r + #(t) fe ete. (6.9) 


the integral-delay representation (6.7) can be developed in terms of successive 
derivatives of x(t). Carrying this expansion to the third order to include the effects 
of radiation reaction one finds the remarkable result 


Bieetecaain? be-2 12. 
Frair(t) = —Mem (#19) & + 5 % (6.9b) 


which agress exactly with the previous equation (2.8), derived on the basis of 
conservation laws.1®) 

The significant new information contained in (6.9b) is the explicit evaluation of 
the electromagnetic mass in terms of the velocity and structural characteristics 
of the electron. Scuorr [30] has carried out extensive investigations of the form 
of (6.9b) for a number of electron models and in particular has verified that in 
the case of the deformable Lorentz electron all the results of the relativistic 
theory are completely reproduced. Repeating this reduction from the appropriate 
modification of (6.7) he showed that 


Mem (x, To) = ees 
: 


in complete agreement with special relativity. Furthermore it can also be shown 
that the coefficient of the % term is precisely equal to (2.10) i. e. the proper relati- 
vistic generalization of mtx. A very important point established by Scnorr in 
this connection is that in the v/c < 1 limit the form of the radiation reaction term 
is independent of the electron model. This is basically the reason that it can be 
recovered at all even in the point electron limit. 

Within the framework of the previous theories it was permissible — in a sense — 
to interpret (6.9b) as being the exact expression for the self-force. Within the 
present development however this appears to be merely the first approximation 
to the structure of a far more complicated theory. It is therefore of particular 
interest to study those features of (6.7) and (6.8) which are not reflected in the 
truncated versions. There are essentially two ways of approaching this problem — 
“One is to remain entirely within the domain of the quasi-stationary approximation 
and to extend the expansion (6.9b) to include the higher terms. The other is to 
go beyond the quasi-stationary approximation altogether and to study the eaact 
consequences of (6.7), for example its behavior in the vicinity of the velocity of 
light and in the “‘super-light”” Cerenkov regime.’’) 


16) More precisely, the quasi-stationary approximation is a necessary but not a sufficient 
condition for deriving (6.9b). The approximation retains its meaning and utility even if 


» > c, but then (6.9b) no longer applies. 
17) Tt is easily seen that (6.9) does not necessarily hold for electrons of energies > 10% ev. 
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We consider first the inclusion of the higher derivative terms: Strictly speaking © 
this involves not only the quasi-stationary approximation but also the additional 
assumption that v<c, and the neglect of all the higher derivative cross terms 
which occur when the expressions (6.9a) are substituted into (6.7). We shall speak 
of these simplifications collectively as the “extended quasi-stationary’’ approxi- 
mation. Since we have already learned from Eliezer’s work (Subsection 6-B above) 
that the inclusion of only a limited number of higher derivatives in (6.9b) is not 
likely to bring about any substantial improvement, we shall immediately extend 
these expansions to include an infinite number of terms: Following the previous 
calculation, the infinite analogues of (6.9a) may be substituted back into (6.7); 
neglecting the products of the higher derivatives, the resulting expressions 
simplify sufficiently so that the remaining integrations can be performed explicitly. 
Formally the result is a differential equation of infinitely high order. The tedious 
details of this computation were actually first carried through by Hererorz [59] 
(see also WiLDERMUTH [4/]) with the final result 


18 e2 eo) 
F seit ae oa s {> Ay or x (6.10) 


where 


Wal = : 6.10 
"= (+2) ln +3] [n + 5]n)) nant 
and ® is the differentiation operator. 
By means of the identity 
~ = exp {hD} —1 (6.11) 


BAe) a [ree ie Le) 


Infinite order differential equations of this kind can always be transformed into 
difference equations. The particular equation (6.10) has however a rather compli- 
cated finite difference form (see subsection F below). It was first shown by 
Pace [15] that this approach leads to much simpler ‘‘natural’’ equations if the 
electron model is taken to be that of a uniformly charged spherical shell. Begin- 
ning with the appropriate modification of (6.7) and again applying the extended 
quasi-stationary approximation one obtains in this case instead of (6.10) the 
expression 


Fe a{> (—2 “Jo «. (6.12) 


2 
315 C \n=1 


Using (6.11), this can now be transformed to 


dee Tea hy 21 : 
Pree = 372 ¢ E ¢ : = |: (6.12a) 


a result which was also later derived by Boum and WEINSTEIN [56]. Comparing 
with our original radiation reaction equation (2.8) we see that now 


Siu 2 
Patz ae E () — 4 (! di: = (6.12b) 
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in contrast to the previous 
PF. a DP e2r ie. 
Rte OS canary 3 


clearly the splitting into inertial and reaction terms has completely fallen away. 
Treating 2rg/¢ as a small quantity, i. e., going back to the partial expansion (6.9a), 
the right hand side of (6.12b) is seen once more to agree with (2.8) in lowest 
order if one chooses 

e2 


2 
m= ———, 
oT 


(6.13) 


The distinction between (6.8) and (2.8) becomes even more marked if one goes 
beyond the quasi-stationary approximation altogether. It was for instance first 
discovered by SOMMERFELD that at ‘‘super-light”’ velocities there is a radiation 
reaction even at constant velocity. By a straightforward modification of the discus- 
sion of Appendix 2 it can be shown that if #(#) = v > c, (6.7) reduces to 


9e2 @ 
Peer = (: (6.13) 


r2 ae 
47% v 


F,.4¢ in this case has the interpretation of being purely a radiation reaction force. 
In fact with the replacement c+ c/n, (6.13) correctly describes the Cerenkov 
losses of a particle traversing a medium of index of refraction n.18) 


D. Finite Difference Theories 


From a logical point of view it is possible to dispense entirely with the scaffolding 
leading from (6.7) to (6.12b) and to adopt the finite difference formulation as the 
fundamental statement of the theory. This point of view and its antecedents has 
been extensively reviewed by CaLprRoua [60]. One can show quite easily for 
instance that it is possible to arrive directly at (6.12b) by merely assuming an 
intrinsic quantization of space and time. In such a lattice manifold, clearly all 
the customary differential expressions have to be replaced by their finite difference 
analogues. In fact, carrying this approach to its logical conclusions, all non- 
constant physical variables must be represented by discontinuous functions. 
CaLpIRoLA himself does not go to quite such extremes but he does present an 
axiom system showing that if one assumes (a) the existence of a fundamental 
interval of time, e. g. 27)/c; (b) the validity of the usual electrodynamics in the 
limit r,-> 0; and finally (c) covariance under Lorentz transformations, then 
(6.12b) emerges naturally as a non-local generalization of the point charge theory. 
_A particulary interesting feature of this approach is that it motivates directly the 
covariant generalization of (6.12b), viz. 


ee fs Wier 
— =~ 44, (s — 8s — 4, (8) xj7(8 — So) x; (s)p = Fy- (6.14) 
saz jie I+ a x (8) i ( 0) i] 

18) Sommerfeld’s formula (6.13) dropped into obscurity because of its apparently anti-relati- 
vistic content. It was not until 1937 that Jorré (quoted in [67]) pointed out that it had 
correctly anticipated the Cerenkov Effect. The elaboration of this correspondence is 
essentially the basis of the phenomenological theory of radiation reaction. [Section 7] 
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The symbols here have their usual meaning as four-vectors; 59 denotes the invari- 
ant interval of proper time corresponding to 2r,/c. This theory is therefore an 


example of a covariant electrodynamics which automatically includes a “‘funda- — 


mental length’. 

CaLDIROLA shows that the runaway features of the point electron model do not 
occur in this finite difference theory (see also subsection F below) and presents 
other arguments showing that the solutions of these equations can always be 
expected to have a dynamically reasonable aspect. Nevertheless there are two 
broad objections which may be brought against too strong a commitment to this 


point of view. The first is that this formulation still does not give a satisfactory — 


account of the energy transfer from particle to radiation field. In particular, as 
CaLpIRoLA shows, the hyperbolic motion is still a radiationless trajectory in this 
theory; in view of the preceding discussions on this point (see especially the 
remarks following (3.7)) it is difficult to accept this without reservation. It 


should also be emphasized that the mathematical structure of (6.12b) is not _ 


nearly as rich as that of its ‘parent’ (6.7): The Cerenkov features and nearly 
all of the non-stationary cases have been completely lost. Nevertheless (6.12 b) 
makes it particularly clear that the vanishing of the external force does not 
necessarily require the vanishing of %. The new solutions which enter here — the 
so-called radiationless internal modes of the electron — will be discussed further 
in the following sections. 


E. Covariant Form Factor Theories 


With the advent of special relativity, the elaborate electron models of ABRAHAM, 
_Lorrntz, and SoMMERFELD were of course completely buried. In fact this inter- 
ment was so complete that when in later years the divergence difficulties of the 
point electron theory were again felt to be particularly acute, the remedies were 
not sought among the old extended electron models but rather among far more 
formal devices — the covariant form factors and smearing functions. The essential 
idea of this new approach was to hold on to the linearity and Lorentz covariance 
of the classical theory but to blur the field-particle interaction sufficiently so as to 
smear out all singularities.) Theories of this type have been studied by Borr 
[13], Feynman [38], McManus [61], and especially LEHMANN [62]; comprehensive 
summaries have been given by H6nt [63] and Rzewusxt [64]. We shall consider 
here only those aspects which have a particular relevance to the radiation reaction 
question. 

In the Hamiltonian form of electrodynamics, the usual interaction between field 
and particle is given by 


Sint =e f [ 0(x — 2(s)) mm (s) Ax (2) dt ds (6.15) 


where z(s) represents a space-time point on the world line of the electron, para- 
metrized by the proper time s, and u, is the four-velocity. In the linear non-local 
theories this is replaced by the ansatz 


Sine = eff (lx — 2(s)]?) u,(s) A, (x) d4x ds (6.15a) 


Several non-local and non-linear theories have also been discussed by Biocutnzey [65]. 


See further Pars and Ustenseck [66]. The significance of radiation reaction in these 
theories is not yet clear. 
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where ¥ is a sharply peaked function which “smears” out the divergences that 
would otherwise be caused by the 6-function. The equations of motion implied by 
(6.15a) are complicated non-linear integro-differential equations. In the non- 
relativistic limit they may however be reduced to the simple form 


+00 
9» 2 ae aw 
Fee) =— 3 8 4 1; dS ((t — tp) re) cures | (6.16) 


where ¥ is related to § essentially via a Fourrer transform. (See [6Z].) One point 
that is immediately evident is that the structure of (6.16) is much simpler than 
that of either (6.6a) or (6.7). Note in particular that the charge volume is inte- 
grated over only once (compare the |x|? in (6.6a)) so that the functional depen- 
dence of F.i¢ can be expected to be more rudimentary in this theory. We shall see 
in a moment that this implies serious drawbacks for (6.16). 


F. Runaway Solutions 


If the non-local or extended electron theories are to represent any improvement 
over the point-charge theory, then clearly they must not contain any runaway 
solutions. The classical investigation of this point dates back to Hrreiorz [59]. 
(See also Wi~tpERMuUTH [4/].) Consider first the case of a uniformly charged 
spherical electron of radius ry. If there are no external forces then (6.3) becomes 


Peg = Mmech - (6.17) 


Assuming the extended quasi-stationary approximation, it is permissible to insert 
(6.10) for Feig, So that 


ie! 9 "0 nr 
18e2 & Cc 


v eee aes Gromes (KG 7/ 
Mech © roc xo |([n + 2] [n + 3] [n + 5] n!) Saree ( #) 


Clearly « = const. is a solution but it is by no means obvious nor even true that 
it is the only solution. To investigate this point further it is useful to make the 
-ansatz 


cA 
=ke — ; 6.17b 
4 = kh exp | oe | ( ) 
From (6.17a) it then follows that 
A 
W (A) exp i = | =a (6.17¢) 
aro 
where ; 
¥ (A) = Mech ae ad ae P(A) (7.16d) 
20 (i 
and 
@ (A) = [A224 442 4+ 475) +44 414 47%. (6.17e) 


Clearly the existence of any solutions of this kind hinges on whether Y (A) has any 
zeros. We state this important result as the 


Qe 
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Herciorz-Witpermutin Theorem: The function Y (A) defined by (6.17 d and e) 
above possesses infinitely many zeros; 4;,7=1, .... These zeros have the 


property that 


SS = 
Re {A;} e Oh AGS Wee - 0. 


The physical significance of this result is that as long as the “mechanical” Mass 
of the electron is not negative, there will be no runaway solutions. The extended 
electron theory therefore not only can cure the runaway difficulties but also 
explains in a very illuminating way why they arise in the first place as ry > 0. 
Basically these “unexpected”? motions of the electron are controlled by two +E 
options: One occurs in the sign of the time — if we reverse the sign of ¢ this 
corresponds to selecting an ingoing solution of Maxwell’s equations (see [3]), 
i.e., energy is fed into the particle and it accelerates. The other option has just 
been revealed by the Herglotz-Wildermuth analysis: If the mechanical mass 
is negative, then the signs are reversed in Newton’s second law and instead of 
the particle experiencing a damping in the face of an opposing force, it accele- 
rates.2°) The limit 7, > 0 enters since we have the upper bound 


Mmech = Mexp — Mem- 


As ry +0, Mem — 00, then clearly mmecn < 0; and this finally is seen to be the 
root of the runaway difficulties ! 

It is interesting to note that this feature has a close parallel in quantum electro- 
dynamics. Since the quantized theory is almost always discussed in terms of a 
perturbation formalism the correspondence is however not usually evident. 
Van Kampen [68] and Norton and Watson [69] have carried out some exact 
calculations with a ‘“‘model’’ quantum-electrodynamics based on a non-relati- 
vistic harmonically bound electron coupled to the radiation field in the electric- 
dipole approximation. By working with an initially smeared electron, it is pos- 
sible to show that as the charge distribution is permitted to shrink to a point, i. e., 
o(x) + 6(x), runaway continuum eigenstates of the Hamiltonian appear preci- 
sely at the point where the “‘mechanical’’ mass of the electron becomes negative. The 
analogy extends even further than this: If one attempts to “regularize” the 
Hamiltonian by altering the commutation relations to suppress this unphysical 
continuum, acausal effects analogous to the classical-pre-accelerations arise.?4) 
We have already indicated several times our preference for drawing a distinction 
between the covariant form factor theories and the extended electron theories. 
The reason for this is that if one carries through an analysis for (6.16) analogous 


*) As 8. CoLEMAN (unpublished) has pointed out this also implies that the particle kinetic 
energy becomes a negative definite quantity — even in the non-point limit. The runaways 
conserve energy simply because they build up the positive field energy (radiation) at the 
expense of negative kinetic energy. See also [90]. 

This point has also been discussed by StrErnwEDEL [14]. In a certain sense in quantum 
theory the situation is even more acute: Since the interactions between the particles 
and the vacuum fluctuations can never be “switched off”, the runaways will always 
persist. The classical adiabatic prescriptions, e. g. [88], then become nugatory. See [90]. 
Another interesting point of the quantum development is that the characteristic classical 
damping constant mc*/e? does not get itself multiplied by factors of e?/hce and therefore 
still remains the natural scale factor. 
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’ to that of (6.17a) one finds the surprising fact that the form factor expression 
(6.16) will in general contain runaway components. This has been verified in 
detail by Bopp [70] and Srrrywepet [14] [70] for a wide class of otherwise coni- 
pletely successful form factor theories. Similar results have also been obtained 
by Irvine [77]. This unpleasant circumstance can best be understood in terms of 
the Herglotz-Wildermuth theorem: The essential point is that most of these 
covariant form factor theories are actually equivalent to a regularization of 
ordinary electrodynamics by a continuum of negative energy meson fields (LEH- 
MANN [62]). The self energy divergences are therefore automatically compensated 
but the electron also acquires a component of negative non-electromagnetic mass. 
According to the Herglotz-Wildermuth theorem one can_ therefore expect 
runaways — and this is indeed what the computations show.22) 


G. Electrodynamic Collective Modes: 
The Excited States of the Electron 


One of the most interesting features of the extended electron theory is its pre- 
diction of the existence of both self-oscillatory and radiationless states of motion. 
This is essentially due to the fact that the basic integral-delay equation 
Feit [x (t)] = 0 of the extended electron theory has associated with it a non- 
trivial eigenvalue problem in the case x(t) ~e’®’. A very simple way of seeing 
this is to go back to the finite difference expression (6.12a) 


ear te 275 
ext | =a = a 
Numer co “| c ) 


which can be derived from the shell electron equivalent of (6.7) in the extended 
quasi-stationary approximation. Clearly any series of the type 


“= > {4 sin (= ) + B;, cos (= i} (6.18) 


=o 0 0 


will satisfy the homogeneous problem F,.i¢ = 0. The physical meaning of this 
result is that an extended charge structure ought to be capable of carrying out 
sustained self-oscillations even in the complete absence of any external forces. Within 
‘the framework of the finite difference theories this is of course an exact result; 
in fact all the trajectories encompassed by the trigonometric representation 
(6.18) can be shown to be strictly non-radiating modes in this formulation [60]. 
In the Herglotz-Sommerfeld theory circumstances are somewhat more compli- 
cated but the structure is proportionately richer and more interesting. The 
appearance of the electrodynamic collective modes in this theory can be under- 
stood in terms of two distinct eigenvalue problems; one corresponding to the 
existence of radiationless modes and the other to the existence of force-free modes. 
We consider these in turn: : | 

(i) A radiationless mode is a sustained charge oscillation which does not emit any 
electromagnetic radiation. The possibility of the existence of these modes was 
first explicitly recognized by EHRENFEST [72] and later by Scuorr [73] although 


22) See also [47]. It is interesting to note that in this respect the “local” requirements imposed. 
by (6.16) are far more stringent than the asymptotic conditions on the Poynting vector. 


380 THOMAS ERBER 


the beginnings of the idea can actually be traced back to SOMMERFELD [74]. A 
concise modern formulation has been given by S. CoremANn (unpublished). The 
essential point is easily exhibited by transcribing the retarded potential expression 
(6.4) into Fourier transform space. Ignoring singular factors, the result is 


Ag. ~SeOr- (6.19) — 


The radiationless condition is then simply 


A, =0 (6.19a) 


which requires that the system oscillate at one of the zeros of o,, the Fourier 
transform of the charge distribution. (Note that (6.19a) applies only to the space 
“outside” the charge distribution.) 

(ii) A force-free mode is a self-sustained charge oscillation which exists even in the 
absence of any external forces. It follows from the previous work that the fre- 


quencies characterizing the force-free modes are the eigenvalues of the functional — 


equation 
Foe (eo) =O: (6.20) 


Roughly speaking for a charge structure of effective radius ~ 10-3 cm these 
frequencies can all be expected to be of the order of 10?3 sect and higher. 

It is important to realize that these force-free modes are not necessarily radi- 
ationless and conversely that radiationless modes are not necessarily force-free. 
Only if the eigenvalue problems (6.19) and (6.20) share a solution does a rigo- 
rously self-sustaining radiationless mode exist. This for example is the case for 
the rotational oscillations of a spherical shell electron [31] [74]. The general situ- 
ation however may be still more. complicated — this for instance will happen if 
either or both of the equations (6.19) and (6.20) have complex roots. In this event 
there may be damped oscillations which nevertheless formally correspond to 
radiationless modes and/or force-free modes which are associated with a definite 
time of decay.?3) 

A further degree of freedom is available in the adjustment of the mechanical 
mass. As WILDERMUTH [4/] has pointed out, if one insists on the auxiliary condi- 
tion Re {/;} = 0, i.e. no damping of the force free modes, then (6.17d) or its 
equivalent introduces still another eigenvalue problem which fixes the value of 
Mmecn- This feature is of course especially interesting in connection with the 
possibility of a many-particle classical field theory. Finally it should be noted 
that in the entire discussion to this point we have been considering only rigid 
charge motions, i. e. the Poincaré stresses have played an entirely passive role. 
Surprisingly high energies may be stored in these collective modes. On a classical 
basis this may be simply estimated by assigning plausible values to the oscillation 
amplitude. For example in the case of the shell electron it can be shown that the 
energy stored in the circulating field components is of the order of one hundred 
times the electrostatic ground state energy. From the quantum mechanical point 
of view it is of course inviting to remove the amplitude ambiguity by assigning the 
energy via Hexe = hvexe; the only remaining imprecision would then lie in the 
*8) These complications require close attention to mathematical details; linearized treatments, 


e. g. [56], [47] can often go astray. A rigorous analysis for one case is carried through in 
Appendix 2. 
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details of the charge structure. If for the sake definiteness one assumes a charac- 

teristic dimension of the order of 1 f, (6.20) again leads to the classical estimate 

exe ~ 10? Hy. Finally one may use a model independent scale factor, i. e. 

2¢?/3mc, the radiation reaction coefficient in the quasi-stationary limit, to 

epee a length. In this case the total energy of the first excited state turns out 
0 be 


3 
mee i a al = 105.54 Mev. (6.21) 


It has been noted previously [56] [75] [60] [76] [77] that this has a suggestive 
resemblance to the muon rest mass 


M,, ='105.65 + .01 Mev. 


7. Absorber Theories of Radiation: The Phenomenological Theory 
of Radiation Reaction 


In the present section we shall discuss two developments which stand somewhat 
apart from the other theories considered thus far — the absorber theory of 
radiation of WHEELER and FEYNMAN [78] and the phenomenological radiation 
reaction theory of GryzBuRG and Erpman [54]. The fundamental point in both of 
these theories is that somehow the presence of absorbers should be crucial to the 
elementary act of the emission of radiation. ,The Wheeler-Feynman approach 
is “fundamental” in its orientation whereas GryzBuRG and Erpman are frankly 
“applied”? — nevertheless these theories have a good deal of conceptual common 
ground albeit there is a lively contrast in some of the ideas involved. We 
consider first the Wheeler-Feynman approach. 


A. Wheeler-Feynman Electrodynamics 


The basic observation of the absorber theory of radiation is the fact that the mere 
presence of an absorber may produce a radiation reaction on the emitting particle. 
This can be most clearly demonstrated in terms of the two body problem of electro- 
dynamics. Consider for example two objects — one far heavier than the other — 
bound to each other through a purely electromagnetic interaction. It was first 
shown by Sync‘) that such a system would lose energy — the particles spiraling 
inwards towards a collision — even without the factoring of the Lorentz force law 
into an external and a damping term (see (2.15)). The system redistributes its energy 
into the radiation field simply through the action and reaction of its two bound 
charges. It is intuitively obvious however that such a process is slow and ineffi- 
cient; indeed Synaz finds that this system collapses at a rate of one thousands 
slower than one which includes explicitly the radiation reaction forces in its 
equations of motion. ie 
The essential point of the Wheeler-Feynman approach now is that it is possible 
to amplify these reactive losses sufficiently so that they become formally identical 
to the ordinary radiation reaction: one simply assumes that the physical solutions 
of the field equation 
vol) 


24) See [79]. A rigorous discussion of this problem has recently been given by Driver [58]. 
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are half the sum of the retarded and the advanced potentials. Formally this obviates 
the sequence of decompositions (2.15) — (2.15e) and automatically inserts the 
correct radiation reaction expression into the equations of motion. Physically 
therefore one has the picture that e. g. a charged oscillator is damped because it 
radiates both forward and backward in time: The radiation travelling forward in 
time causes electrons in the absorber to move — these in turn emit radiation both 
forward and backward in time. The radiation travelling backward in time reaches 
the oscillator at the same time that the original radiation was emitted and exerts 
a retarding force on the oscillator — this is the radiation reaction. The total 
field of oscillator and absorber combine to give the fully retarded field which corres- 
ponds to ‘experience’. 

Naturally the appearance of the acausal advanced potentials is the most startling 
feature of this theory, and WHEELER and FEYNMAN have gone to great lengths 
to show that in the case of complete absorption the results are formally identical to 
those of the conventional radiation reaction theory. (See also [9/]). In particular all 
of the details of the absorber structure are shown to disappear from the final 
formulas; a result which in view of the discussion of the next section is perhaps not 
entirely fortunate. Moreover it is no easier in this formulation than in any other 
to see in detail how the radiation reaction of certain many body problems e. g. 
a ring charge drops to zero as one makes the transition to the continuum limit. 


B. The Theory of GryzBure and ErpMAn 


The phenomenological radiation reaction theory of GrnzBuRG and ErpMAN [54] 
[80] is an attempt to give a unified description of the dynamical consequences of 
certain interactions of charged particles with matter. It will be helpful in fixing 
the ideas involved if we first give a brief qualitative summary of some represen- 
tative experimental situations. It should be noted that Bremsstrahlung and other 
microscopic phenomena are ignored in the subsequent discussion. 

(i) The ordinary Cerenkov effect is one of the most elementary examples of a 
departure from the 


PRS Ne 
Preact = rs aor (7.0) 


law of the point-charge vacuum theory. The correct expression for the radiation 
reaction in this case is 


ee. C 2 
Erne = — et eft — (azar) | oda. (7.1) 


v>eln 


If the index of refraction n is a sufficiently slowly varying function of w, and if the 
integral is cut off at w = ¢/r9, it can be shown that (7.1) reduces to the expression 
previously given by SoMMERFELD [see (6.13)]. It is interesting to note that such 
a formula may continue to apply even if the particle itself is moving in a vacuum: 
This for instance will be the case if the motion occurs in a narrow channel embedded 
in a medium with index of refraction n > 1. A futher extension of this description 
is useful for motions which occur in a waveguide or near some other periodic 


structure ; for example, in a traveling wave tube used for the generation of milli- 
meter microwaves. 
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~ (ii) The complex Doppler effect (Cerenkov self-excitation): Consider a system 
possessing a natural radiating frequency w, (measured in the rest system) moving 
with a velocity v through a medium characterized by an index of refraction n (a). 
The Doppler shifted frequency wg observed at an angle @ with respect to the 
direction of motion is then in general determined by the equation 


Se al te Py 
‘ |1 — Bn(a) cos O|’ ies eet (7.2) 


Suppose now that the auxiliary inequality 
v cos O > w(a) (7.3) 


where w denotes the group velocity, is satisfied for a value of «wq; it then follows 
automatically that for a fixed © there will actually be at least two Doppler 
components, 1. e., @ga—> mq and wy which are solutions of (7.2). This con- 
dition is the Complex Doppler Effect.?°) 


Suppose furthermore that the system velocity v is sufficiently great so that 
Bn(w) > 1. (7.4) 


In this case, in addition to the ‘‘doubled”’ range of Doppler frequencies the system 
will also be able to emit Cerenkov radiation. Let the opening angle of the Ceren- 
kov cone corresponding to wg be Qg, i.e. 


Og = sec [Bn (wa)]. (7D) 


Then for all Doppler frequencies emitted into the cone 0 < 0, we will have the 
remarkable circumstance that the emission of radiation will be accompanied by 
a simultaneous excitation of the system. (Cerenkov Self-Excitation!) Although 
such a result is at first sight surprising it is in fact a simple consequence of the con- 
servation of energy and momentum for the emission of radiation from a quasi- 
bound (i.e. n > 1) system. Quantitatively this is described by 


AU = —E(a) [1 — Bn(a) cos O] (7.6) 


where E(w) is the Cerenkov energy at frequence m and AU is the corresponding 
increment of internal energy. Clearly if Bn(w) cos O > 1, then AU > 0; and this 
is the condition for Cerenkov self-excitation. 

This self-excitation has the following dynamic consequences: Suppose that a 
harmonic oscillator of frequency w, vibrates along the z-axis and simultaneously 
is in translation such that vay = v, > c/[n(@) cos @]. In this event the emission 
of radiation to the regions outside the Cerenkov cone will be associated with a 
damping of the oscillator vibrations; on the other hand the radiation into the Ce- 
renkov cone will be associated with an increase of the amplitude of oscillation. 
There is of course no question of a violation of the conservation of energy here 
since both the oscillation and the radiation occur at the expense of the kinetic 
energy of translation. Nevertheless this presents the novel feature of a radiation 


25) I M. Frank [87]. Necessary conditions for the occurrence of the Complex Doppler Effect 
have been derived by BarsuKov [82]. 
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reaction force formally decreasing the damping of a radiating oscillator. As empha- 
sized by GrnzBure this circumstance is of practical importance in considering the 
beam stability of “super-light”’ bunches in a plasma accelerator.*) ‘aid 

(iii) The Inverse Complex Doppler and Cerenkov Effects [84]: Radiation 
reaction becomes an even more complicated affair if one admits the possibility 
of negative group velocities. This is a matter of some practical importance in 
media having anisotropic and gyrotropic properties as well as in regions of 
high anomalous dispersion. A simple illustration is provided by an idealized me- 
dium with «<0 and mw <0. Since n = (eu)'/* remains a positive quantity in 
such a medium undamped electromagnetic waves may still be propagated. How- _ 
ever phase and group velocities will now have opposite directions since the plane » 
wave condition requires that 


S=C|EPk (7.7) 


and for this case ¢ = —|e/u|'. (S, k, and # have their usual significance as 
Poynting vector, wave vector, and electric field respectively.) 

An interesting peculiarity that makes its appearance here is the necessity for 
using advanced potentials in the construction of solutions of Maxwell’s equations. 
To preserve the correct sense of the energy flow the sign of ¢ must be reversed in 
tandem with the extra minus sign acquired by k& because of (7.7). Physically this 
means that the conventional picture of the Cerenkov cone is reflected into its 
mirror image, i.e. the Cerenkov cone makes an obtuse angle with the particle 
motion. A particular consequence of this is that the transition radiation emitted. 
in the passage of fast particles from vacuum to medium will in this case go back- 
wards into the vacuum. If now in addition there are Doppler frequencies present 
due to internal degress of freedom, then the low frequency components will be 
radiated forwards while the high frequency components will go backwards. (In- 
verse Doppler Effect!) Finally if we have all three conditions present simulta- 
neously, i. e. Cerenkov effect, Doppler emission, and negative group velocities; 
then it may happen that the anomalous Doppler components will outweigh the 
normal components, so that the total emission will in fact be associated with a net 
excitation of the oscillation. The radiation reaction work in this event is a purely 
positive quantity. It should be noted that this is in sharp contrast to the assump- 
tions of Section 5. 

(iv) An idealized isotropic plasma has a refractive index 


|" 
i= |e 


mw? 


for transverse electromagnetic waves. An oscillator of natural frequency w2<4e2N/m 
embedded in this plasma will therefore be unable to lose energy by radiation. 
This brief summary should make it clear that the construction of a phenomenolo- 
gical radiation reaction theory endowed with sufficient mathematical flexibility 
poses a formidable technical problem. The most natural idea would be of course 
to start from some fundamental “vacuum” theory —the Freact = (2€2 /3c3) & version 
for example — and then to work forwards towards a macroscopic description via 
suitable averaging processes. It should be realized however that this cannot 
simply be a matter of imitating the procedures of ordinary electrodynamics: The 


*6) If wo = 0, (7.3) still is relevant as a criterion for the onset of the Complex Cerenkov Effect, 
AGRANOVICH et al [83]. 


| 


The Classical Theories of Radiation Reaction 385 


usual transition (Microscopic Maxwell Equations) = (Phenomenological Max- 
well Equations) merely involves the replacements 1 => je], 1= ||u|| in the 
constitutive relations; the essential structure of the differential equations remains 
intact. In the radiation reaction case it should however be obvious, e. g. from the 
contrast between (7.0) and (7.1), that more drastic alterations will be required.?7) 


GrxzBurG and Erpman [54] meet this problem by going back to the Sommer- 
feld extended electron picture. The radiation reaction is once again visualized 
as originating in a net force resulting from the action of one part of the electron 
on another. It has already been shown (Section 6-C) that this approach gives rise 
to a mathematical structure of considerable complexity. The crucial departure in 
the Ginzburg-Eidman development is that the information on the macroscopic 
structure of the medium is inserted into the formalism by assuming the validity 
of the constitutive equations 


D=|el|H, B=|u| Hz 


an the interior of the electron. This leads precisely to the previous Fourier self-force 
representation (6.6a) except that now the index of refraction appears explicitly 
in the integrand. In a detailed discussion GrnzBuRG and ErpmAn show that this 
recipe is indeed capable of reproducing some of the behavior expected in the 
situations enumerated above. The spirit of the treatment is however completely 
pre-relativistic — an especially unpleasant feature considering the high velocities 
encountered in the Cerenkov effect. GinzBurG and E1pMaAN attempt to minimize 
this embarrassment by making the transition @(r) + 6(r) as soon as feasible; the 
subsequent divergence is then absorbed as a mass renormalization. But this pro- 


cedure is itself of course neither covariant nor model-independent. A counter 


example was provided long ago by SOMMERFELD [3] in his proof that the spherical 
shell electron would suffer infinite radiation reaction under Cerenkov conditions. 
It may be concluded that the formulation of a completely satisfactory macroscopic 
radiation reaction theory still remains as an unsolved problem. 
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Appendix 1! 
Inequalities Connected With the Radiation Equation 


We consider the modified Fowler-Emden equation 


ee U0 tee Ai.) 
j 5 ( 


27) Another important point that arises here is the “Synge dissipation” [79]. If indeed two 
charges can emit radiation even in the absence of any explicit damping terms, then one 
must anticipate the possibility of cooperative dissipation mechanisms in a macroscopic 
theory. (I am indebted to Prof. Prigogine for an enlightening discussion on this point). 
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on the closed interval Q.: 0<t<t,, under the assumptions 


(i) g(t) continuous on 2, 
(ii) W(t) increasing and non-negative on Qo, 
(iii) y(t) and y(t) non-negative on Q¢. 


and with the initial conditions 


(iv) y(0) = 0; y(0) =a, >0. 


From (i) — (iii) it follows that y (é) is a concave function. In particular from (iv) we 


conclude that 
y(t) < aol. 
Inserting this into (A1.1) we obtain the inequality 
a Ps apt 


Integrating from 0 to ¢, there results 


t 
1 W 
— (i) + a% > eae dt; 
0 


where the existence of the integral is guaranteed by (A1.3) and (i). 


But now 
t t 
W W W Ww 
ee ares 
0 0 


which again is guaranteed by (i) and the mean value theorem. 
Therefore 


WS 
Since a > a) — y(t), this yields the important result 
Wey nae 
We may now perform another quadrature: From (A1.6) it follows that 


ULCSY. cae 
Ao 


v= [ (Ta 


and (A1.7) assures the existence of the integral. 


where 


(A1.2) 
(A1.3) 


(A 1.4) 


(A .15) 


(A1.6) 


(A1.7) 


(A1.8) 


(A 1.9) 


The crucial point of the argument now becomes apparent: If the original inequality 
(A1.2) is taken through this cycle, then the end result (A1.8) is an inequality 


a 
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which is an improvement over the input. The expectation therefore is that the 
entire process can be iterated with successive improvements. Let us consider the 
next cycle: 


Putting y,(t) = ayt — W/ay; then from (Ai.1) it follows that 


W 
= nih (A1.10) 
Integrating, i 
t 
Ww 
—yH +a Hee dt; (A1.11) 
© WY 


where now 
t t 


{( di = — +{ (3) y, dt > ua (A1.12) 
; Yy Yy , Yi yy 


and we note that the boundedness of the derivatives and the mean value theorem 
guarantee stepwise the existence of each subsequent parts integration. 
Therefore 


W(t) < ayy, (0): (A113) 
and integrating again 
€ 
ayt — y(t) >| (=) dt. (A1.14) 
‘ Yn 
But now : 
gees (1 : at (A115) 
Yy At aot 
and 
t t 
7 2 172 72 
Ue Se 1 : 1 WwW it [iv (A1.16) 
(2 2 2% a a ae od) 
Therefore the final result of the second iteration is 
vo 
Ae aE age aaa JNA Ty 
OES ede (A1.17) 


which, as expected, is a further improvement over (A1.2). It is now clear that 
this pattern may be repeated n-fold. On the n' iteration one obtains 


n Ww 
y(t) < aot { sei Da Bs 


k=1 
or, proceeding to the limit, 


y (t) < aot e — exp pal (A 1.18) 


This is the result stated in the text as (5.12). Since ¥ > 0 for ¢ > 0, (A1.18) may 
also be considered a transcendental inequality for W ; the solution is elementary and 
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yields (5.13b) of the text. The n® analogues of (A1.7) and (A 1.13) are the W (ty 
restrictions (5.13a). It is interesting to note that this process yields simultaneous 


restrictions on y and W. 


Appendix 2 
Sommerfeld’s Self-Force Representation 
Theorem: The trajectory defined by a(t + 2ry/c) = x(t), |x| <e, for a uniformly 
charged electron of radius 79, is neither radiationless nor force-free. 
Proof: Let 
e(t, t) er AY where 7 = 2) — a(t =a) (A2.1) 
Then 
2r 
e(t, %) —2% =0; m=: (A2.1a) 


This solution is unique since de/0t > 0 and e(t, 0) = 0. 
Similarly if 

é (t,t) = ct — 7, 
then i 
(t,t) > 0 and F(t, 7,) —2r,=0. (A2.1b) 


This establishes that the equations ct + T = 27); |ct —T| = 2r, both have 
the unique root t = Ty). ; 
The self-force representation (6.7) may therefore be rewritten in the form 


3 2 
Foot = — 55 (Li + In + Is} (A2.2) 

where ; 
he — fac a (f(er + T) —f(et —T)} (A2.2a) 

0 
° @e 
Tees fa + ag her + 1) — fer — 7)} (A2.2b) 
0 

fos ae yaks 

I= f dr 0) Z Her $9) + er — 1). (A2.20) 
0 


Inserting the explicit form of f [see (6.7a)] we have 


To 


4 
| <e fax {7 ‘Sr — 6(cr)* + Sroer| + 


0 


ct \2 “ fs 1 
+9|(=) -2| ra cor | (A2.3a) 


| 
/ 
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To 


i= c fac {7 ee — 6(ct)? + 87, (en| + 
5 0 


cr \? 1 
+9[3(Z)—s]40+[ a reat 
= fa i+) {l? ce = Sines sr 4 q alli (A2.3¢) 
if 0 0 


The {-1 terms give rise to a divergence which cancels out when the integrals 
are combined. It is now convenient to transform J, by a parts integration; this 


results in 
(er)? | 2 
[= = 3 
; c [ax |z|° 2a an }. (A2.3d) 
i 
Collecting all the terms we finally obtain 
I, I, + 1, = 4e " i ieee — 73 we A2.4 
“al tide aon : a Gory ee 
: 0 


(Since  ~ O (vr,/c), a lear approximation would now be tantamount to drop- 
ping the {? term.) 

From (A2.4) we read off that Fee (¢ + t)) = Fseig (t); % is of course also periodic 
with period 7». It is therefore sufficient to consider the work done over a period: 


W (to) = f Feerr(t) &(¢) dt. (A2.5) 
0 
In carrying out this computation one encounters three different types of integrals: 
J, = f dt {dre (A2.6a) 
0.8 | 
Jn = [dt {dravT (A2.6b) 
Ont 7G 
Jn = f dt {drag (A2.6c) 
0 0 


With the help of the auxiliary integrals 


a [drat —2 =0 
0 


aren x(t) =0, n>0 (A2.7) 
0 
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i 


it is easy to show that 
Sia re ge (A2.8) 


However J, has a non-vanishing component. Let 
Jp = IP + IP =f dt f dr[PeQeO) — fat fdr [x2 w(t —r) &(é)]. (A2.9) 
0 0 0 0 


Then J® = 0 because of (A2.7). J? may be transformed by successive integ- 
rations by parts; first with respect to ¢, then with respect to Tt: 


IP = —t4 [ dta() + IP; 
0 


To To 


J@ =2[ dt { dr [re(@)x(t—z)]. (A2.9a) 
0 0 


Now we note that W (z,) is invariant under the transformation x(#) > x(t) +- 
const. There will therefore be no loss of generality in assuming that y = x also 
defines a periodic function. JS may then again be integrated by parts just as 
in the transition from (A2.9) to (A2.9a). 


To 


J® = ro f aty(t) y (t) — far fartyoat —t)|=0. (A2.9b) 
0 0 


0 


Finally therefore we obtain 


3 2 
W(t») = ——  f aearey 20, (A2.10) 
2r5 
A fortiori this proves that 
Peg == 0. (A2.11) 
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